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ABSTRACT 


In  earlier  portions  of  the  present  research  program,  a  precise 
mathematical  model  was  developed  for  the  estimation  of  sound 
absorption  over  a  ground  cover  layer.  The  usefulness  and 
validity  of  this  analysis  was  confirmed  by  r-“sults  obtained 
through  a  laboratory-scaled  experiment.  In  the  present  study, 
the  above  mentioned  mathematical  model  is  applied  to  the 
prediction  of  sound  absorption  over  various  types  of  ground 
covers  composed  of  natural  vegetations.  Through  an  extensive 
review  of  literature  in  agriculture  and  forestry,  the  physical 
structure  and  relevant  mechcniccl  properties  of  ground  cover 
canopies  have  beer,  determined.  The  acoustical  properties  of 
the  ground  cover  itself  has  also  been  estimated.  Based  on  these 
results,  ground  absorption  spectra  over  various  assumed  ground 
cover  conditions  have  been  calculated.  The  trends  of  the 
prediction  have  been  compared  with  experimental  results  and 
the  limitations  of  the  present  study  are  discussed. 
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1 .0  INTRODUCTION 


In  earlier  parts  of  the  present  study  of  sound  absorption  by  natural  ground  cover,  a 
mathematical  model  has  been  developed  for  estimating  the  sound  absorption 
spectrum.  The  usefulness  and  validity  of  this  model  was  confirmed  by  results  from  a 
laboratory-scaled  experiment  of  sound  attenuation  by  simulated  ground  cover 
(Reference  1 ).  In  this  model  the  ground  cover  is  represented  as  a  layer  of  acoustical 
material  of  finite  thickness.  The  concept  is  an  extension  of  a  representation  chosen 
by  Ingard  (Reference  2),  where  the  ground  or  ground  cover  is  des  jnated  as  a  semi- 
infinite  medium.  The  analysis  is  straightforward,  and  should  be  very  accurate  in 
the  far-field. 

The  peak  absorption  frequency  and  the  general  shape  of  the  sound  absorption  spectrum 
depend  upon  several  principal  parameters:  the  height  of  the  sound  source  above  the 
ground;  the  thickness  of  the  ground  cover;  acoustical  properties  of  the  ground  cover; 
and  the  normal  acoustical  impedance  on  the  surface  of  the  semi-infinite  ground.  In 
this  report,  a  detailed  study  of  the  mechanical  structures  of  natural  ground  cover 
such  as  forests  and  field  crops,  has  been  undertaken.  According  to  their  mechanical 
properties,  the  acoustical  properties  of  natural  ground  cover  are  estimated.  Thus,  it 
is  feasible  to  make  realistic  predictions  of  sound  absorption  spectra  of  natural  ground 
covers  over  a  considerable  range  of  conditions.  The  results  are  presented  in  both 
graphical  and  tabulated  forms.  tAfith  increased  knowledge  of  the  acoustical  properties 
of  the  ground  cover,  refinements  of  the  predicted  values  can  be  made  in  the  future. 
However,  the  present  results  are  considered  to  be  suitable  for  preliminary  estimates 
of  ground  absorption  of  sound  in  field  conditions. 

In  the  remainder  of  this  report,  a  condensed  description  of  the  analysis  and  the  resulting 
equations  are  given  in  Section  2.0;  a  general  description  of  the  structure  cf  different 
types  of  forest  and  field  crops  is  presented  in  Section  3.0;  the  detailed  estimation  of 
acoustical  properties  of  natural  ground  covers  is  discussed  in  Section  4.0;  the  pre¬ 
dictions  of  ground  absorption  spectra  under  a  range  of  given  conditions  are  presented 
in  Section  5.0;  and  finally,  the  limitations  of  the  present  study  and  other  general 
discussions  are  included  in  Section  6.0.  In  addition,  a  FORTRAN  computer  program 
which  performs  the  calculation  of  the  sound  absorption  spectrum  is  documented  and 
listed  in  the  Appendix. 
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2.0  MATHEMATICAL  ANALYSIS 


The  precise  analytic  nature  of  sound  attenuation  near  a  boundary  with  known  acoustic 
impedance  was  first  made  clear  in  a  series  of  studies  by  Rudnick  (Reference  3),  and 
Ingard  (Reference  2).  The  main  application  has  been  the  estimation  of  ground  attenu¬ 
ation  effects  on  sound  propagation  in  the  atmosphere.  Some  subsequent  analytical 
studies  and  experiments  have  further  explored  the  details  of  this  phenomenon,  in  these 
studies,  the  acoustic  media  above  and  below  the  boundary  plane  are  assumed  to  be 
semi-infinite.  Constant  values  are  ascribed  to  either  the  impedances  of  the  two  media 
or  the  normal  impedance  of  the  boundary  itself.  However,  in  many  situations  with 
practical  importance,  the  boundary  between  the  upper  and  the  lower  semi-infinite 
acoustic  media  is  not  a  simple  plane,  but  a  porous  layer  with  finite  thickness.  It  is 
natural,  then,  to  investigate  the  wave  attenuation  characteristics  near  such  a  com¬ 
posite  boundary. 

Owing  to  the  special  nature  of  wave  reflections  at  near  glancing  angles  of  incidence, 
it  is  not  possible  to  assign  a  constant  value  of  impedance  to  the  entire  composite 
boundary.  A  layered  media  representation  becomes  necessary  in  this  case.  It  is 
necessary  to  specify  both  the  acoustic  impedance  and  the  wave  transmission  constant  of 
the  porous  transition  layer.  A  new  analysis  is  therefore  required  for  estimating  the 
magnitude  and  characteristics  of  the  ground  attenuation  of  such  a  layered  boundary. 

In  addition  to  the  analysis,  an  experimental  study  was  undertaken  in  the 
previous  year  (Reference  1)  to  determine  the  attenuation  characteristics  of  a  layered 
boundary.  The  results  of  the  experiment  were  intended  for  the  verification  of  the  theory 
as  well  as  for  obtaining  a  separate  view  of  the  problem  from  an  independent  approach. 
The  experiment  was  performed  under  laboratory  conditions  so  thar  the  ground  attenu¬ 
ation  effects  were  studied  without  the  uncertainty  of  other  complications,  such  as  wind 
refraction  and  turbulent  scattering .  Overall,  the  study  found  significant  departure  in 
several  aspects  of  the  attenuation  characteristics  of  a  layered  media  from  those  of  a 
simple  boundary.  Thus,  one  may  find  the  results  useful  in  dealing  with  a  variety  of 
practical  problems  where  a  layered  representation  of  the  boundaries  is  warranted. 

Geometrically,  the  space  is  assumed  to  be  divided  into  three  layers.  The  top  semi¬ 
infinite  layer  is  assumed  to  be  air,  which  has  a  density  of  pQ  ,  speed  of  sound  cQ  , 

and  acoustic  impedance  of  p^  c  .  The  middle  layer  is  assumed  to  be  a  porous 

material.  Its  density  and  speed  of  sound  are,  in  general,  complex  quantities.  In  other 
words,  it  has  a  complex  acoustic  impedance  such  that  a  plane  wave  transmitted  from 
the  air  into  this  layer  will  be  refracted  into  the  layer  with  a  phase  shift,  and  will  be 
attenuated  as  it  propagates  through  this  material.  A  third  medium,  which  represents 
the  ground,  occupies  the  lower  half-space.  To  simplify  the  analysis,  a  constant 
normal  impedance  is  prescribed  at  the  interface  between  the  middle  layer  and  the 
semi-infinite  ground. 

\ 


2 


As  pointed  out  in  previous  studies,  the  simple  ray  acoustics  approach  cannot  account 
for  the  observed  ground  attenuation  phonomenon,  and  a  more  rigorous  mathematical 
analysis  must  be  followed.  In  the  present  study,  the  approach  of  Ingard  (Reference  2), 
together  with  coordinate  systems  and  symbols  in  that  paper,  is  adapted. 

A  spherical  wavefront  which  originates  from  a  point  source  can  be  represented  as  an 
integral  of  its  plane-wave  elements 
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In  the  above  integral,  the  vector  (k  ,  k  ,  k  )  inderfifies  the  wavenumber  of  a  plane- 
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wave  element,  and  h  denotes  the  height  of  the  sound  source  above  the  top  of  the 
layered  boundary.  The  reflection  of  the  primary  wave  at  a  boundary  can  then 
be  represented  as 


2 n  n/2J>  I  ao 

pr  -  (£)  f  »f  .:V+V+S<h+z)]R<(M>sinede  (2) 

o  o 

where  R (0, 8)  is  the  plane-wave  reflection  coefficient.  The  reflection  coefficient 
R(0,0  )  is  a  function  of  0  and  6  . 

It  is  more  convenient  for  the  purpose  of  integration  to  write  Equation  (2)  in  a  new 
spherical  coordinate  system  where  the  principal  direction  of  the  reflected  ray  is 
chosen  as  the  reference  axis  (Reference  2).  The  new  angular  variables  are  defined 
as  0  ana  q  •  Equation  (2)  can  now  be  represented  as 

r2n  rn/  2+iao 

Pr  =  {~5tt  )  J  J  e'kr2COSn  R(0,q)sinqdq  (3) 

(f  0 

An  integral  of  this  form  can  be  evaluated  by  using  the  method  of  steepest  descent 
(Reference  9)  in  the  acoustic  farfield,  where  the  value  of  kr^  is  large  compared  to 

unity.  Along  the  paths  of  steepest  descent  in  the  complex  q  plane,  a  new  variable 
can  be  defined  such  that 


cos  q  =  1  +  it 


(4) 
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where  t  is  real  and  positive.  The  reflected  wave  can  be  then  written  as 
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which  has  the  form  of  a  wave  originating  from  an  "image  source"  located  at  o  distance 
h  below  the  top  of  the  layered  boundary,  with  a  variable  strength  Q. 

For  an  arbitrarily  given  function  of  R(t£,8),  Equation  (5)  can  only  be  integrated 
approximately.  A  first  asymptotic  approximation  can  be  given  as 
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where 


y  =  cos  9 
0  0 


R‘  and  R"  are  the  first  and  second  derivates  of  R  with  respect  to  y  =  cos  0  .  Equation 
(6)  serves  as  a  starting  point  for  the  present  analytical  investigation  into  ground 
attenuation  due  to  layered  media. 

It  remains  here  to  determine  the  plane-wave  refraction  coefficient  R (0)  for  a  boundary 
with  a  layered  configuration.  The  over-all  reflection  coefficient  accounts  for  the 
wave  reflection  at  the  top  of  the  middle  layer,  as  well  as  the  wave  that  is  transmitted 
into  the  middle  layer,  relfected  by  the  ground,  and  returned  into  the  air.  Hence, 
the  over-all  reflection  coefficient  of  the  layered  boundary  can  be  given  as 


R(0)  = 
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where  5^  is  the  specific  admittance  rotio  of  the  middle  layer  with  respect  to  air, 

B  is  the  specific  admittance  ratio  at  the  interface  of  the  middle  layer  and  the  ground, 

2 

and  0  is  the  refraction  angle  in  the  middle  layer.  The  function  0  accounts  for  the 

phase  difference  between  the  two  reflected  wave  components.  This  phase  shift  is 
caused  by  path  difference  and  the  wave  transmission  characteristics  of  the  middle 
layer. 

The  insrantaneous  value  of  sound  pressure  in  the  farfield  can  now  be  determined  as 


(i- y  2)R,,(y  )  -  r  R ■  (r  ) 
o  0  0  c 


The  derivatives  of  R(0)  with  respect  to  cos0  can  be  obtained  from  Equation  (7).  By 
substituting  the  known  expressions  for  R(0),  R 1 ,  and  R"  into  Equation  (6),  an  explicit 
expression  can  be  obtained  for  the  reflected  wave.  The  attenuation  of  sound  near  the 
layered  boundary  can  now  be  obtained  by  simply  adding  the  incident  and  the  reflected 
sound-pressure  fields. 

The  algebraic  expressions  involved  in  computing  the  sound-pressure  field  in 
the  upper  half  space  are  straightforward  but  relatively  bulky.  Therefore,  the  results 
have  been  programmed  for  computer  calculation.  The  asymptotic  approximate  solution 
is  very  accurate  for  computing  sound-pressure  levels  in  the  farfield,  i.  e.,  points 
that  are  more  than  a  few  wavelengths  away  from  the  sound  source.  For  the  study 
of  ground  attenuation  effects,  this  is  an  insignificant  restriction.  In  the  computing 
program,  all  of  the  geometrical  and  acoustical  parameters  can  be  varied  independently. 

In  particular,  the  specific  admittance  ratios  B  ,  B  and  the  refraction  index  n  are 
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assumed  to  be  complex  numbers . 


3.0  THE  STRUCTURE  OF  FORESTS  AND  FIELDS 


The  structure  of  natural  ground  cover  is  determined  by  a  wide  variety  of  factors. 
Climate,  terrain,  soil  condition,  rainfall,  plant  species,  ond  the  plant  community 
regeneration  process  all  come  into  effect.  For  the  present  study,  it  is  important 
to  know  the  values  of  parameters  such  as  height,  foliage  density  in  the  canopy, 
stratification  of  the  canopy,  nature  of  the  undergrowth,  and  seasonal  variations. 

Many  of  these  properties  have  been  studied  extensively  in  forestry  literature. 

However,  since  each  forest  being  studied  has  its  own  special  features, descriptions 
are  very  difficult  to  make  (References  4  through  10). 

In  this  section,  the  general  appearance  and  common  variation  of  ecological  conditions 
of  three  major  categories  of  ground  cover  vegetations  will  be  discussed.  These  main 
categories  ce: 

•  The  temperate  forest 

•  The  tropical  rain  forest 

•  The  grain  crop  and  the  grass  fields 

The  outward  appearance  and  the  fine  structure  of  ground  cover  are  entirely  different 
for  each  of  these  broad  categories.  Further  classification  with  each  category  will  be 
defined  in  the  course  of  the  discussions. 

3.1  The  Temperate  Forest 


In  the  temperate  forest,  the  plant  community  is  commonly  dominated  by  individual 
trees  from  only  a  few  species.  A  single  temperate  forest  rarely  contains  more  than 
twenty  species  of  plant  life  forms.  Most  of  the  time,  the  ranking  of  the  plant 
species  and  their  roles  in  the  forest  are  clearly  recognizable.  The  trees  in  a  tem¬ 
perate  forest  can  be  classified  into  five  groups:  standard  trees;  dominant  trees; 
codominant  trees;  dominated  trees;  and  suppressed  trees.  Their  descriptions  are 
given  below. 

Standard  Trees:  These  are  isolated  large  trees  in  the  stand.  Their  trunk 
diameters  are  significantly  larger  than  average  and  the  crowns  are  fully 
developed. 

Dominant  Trees:  As  a  rule,  these  trees  form  the  main  part  of  the  stand, 
and  have  relatively  well  developed  crowns. 

Codominant  Trees:  These  trees  have  fairly  normal,  but  comparatively 
weak  and  narrow  developed  crowns. 

Dominated  Trees:  The  crowns  are  more  or  less  stunted  with  one  side 
developed,  or  suffer  severe  pressure  on  one  or  more  sides. 
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(a)  Trees  with  crowns  in  the  middle  story,  their  heads  mainly 
free  but  in  most  cases  completely  surrounded. 

(b)  Trees  with  crowns  partly  in  the  undergrowth. 

Suppressed  Trees:  These  trees  are 

(a)  Trees  with  crowns  capable  of  survival. 

(b)  Trees  with  crowns  dying  or  completely  dead. 

The  main  foliage  canopy  is  formed  by  the  crowns  of  the  standard  trees,  dominant  trees, 
and  the  codiminant  trees.  Trees  in  these  three  classes  have  approximately  the  same 
height;  their  differences  are  mainly  their  sizes  and  the  crown  development.  The 
bottom  of  the  canopy  is  clearly  defined  in  most  cases.  Undergrowth  in  a  temperate 
forest  is  usually  sparse  and  orderly,  and  is  composed  of  mainly  seedling  trees,  sapling 
trees,  other  small  trees,  and  shrubs.  The  height  of  small  trees  and  shrubs  is  somewhat 
lower  than  the  bottom  of  the  main  canopy.  In  most  of  the  temperate  forest,  there  is 
only  one  main  canopy.  Only  on  rare  occasions  are  there  two  stories  of  canopy. 

On  the  floor  of  the  forest,  there  is  commonly  a  shallow  layer  of  fallen  leaves  and  twigs, 
together  with  some  growth  of  herbaceous  plants.  Otherwise,  there  should  be  no  major 
obstacles  to  interfere  with  clear  passage  on  foot.  Visibility  through  the  tree  stands  is 
variable,  depending  mainly  on  the  height  of  the  lowest  level  of  branches  on  a  tree, 
and  the  seasonal  manifestation  of  foliage  and  floration  in  the  forest  trees  and  the 
undergrowth. 

There  are  two  main  types  of  forests  in  the  temperate  zone:  deciduous  and  evergreen. 
Most  of  the  deciduous  trees  are  broad-leaf  trees.  However,  a  few  species  of  needle- 
leaf  trees,  such  as  eastern  larch,  are  also  deciduous.  In  a  deciduous  forest,  seasonal 
variation  of  foliage  color  and  density  can  provide  very  striking  changes  of  the 
appearance  of  the  forest  stand.  A  mature  deciduous  forest  reaches  a  height  of  approxi¬ 
mately  80  ft.  The  bottom  of  the  main  canopy  is  30  to  40  feet  above  the  ground,  and  the 
closed  canopy  is  about  30  to  50  feet  in  thickness.  A  typical  deciduous  tree  model  is 
shown  in  Figure  1  (b).  Owing  to  the  variation  of  weather  and  water  supply,  the  foliage 
density  in  the  main  canopy  can  change  not  only  with  the  season,  but  also  from  year  to 
year. 

The  dominant  species  for  forest  stands  in  the  same  general  geographical  region  can  be 
different  from  location  to  location.  The  overall  pattern  may  appear  to  be  a  mosaic  of 
splendid  proportion.  It  is  interesting  to  note  that  the  seedlings  and  saplings  in  a 
temperate  forest  stand  may  not  belong  to  the  same  species  as  ihe  dominant  or  codominant 
trees  in  the  same  stand.  This  is  because  the  saplings  of  the  dominant  species  may  not 
be  able  to  survive  the  shade  under  the  forest  canopy.  Hence,  young  trees  of  a 
persistent  and  shade  tolerant  species  will  flourish,  and  eventually  take  over  as  the 
dominant  species  in  the  next  succession.  Such  a  sequence  of  change  in  the  dominant 
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Figure  1  .  Model  Trees  in  Temperafe  Foresfs.  All  dimensions  are  in  feet. 
(Reference  8) 


species  has  long  been  observed.  For  example,  red  maple  will  succeed  blackgum,  and 
be  replaced  by  beech  (Reference  9). 

In  a  temperate  evergreen  forest,  the  dominant  species  are  needle-leaf  trees.  These 
species  are  generally  called  conifers.  A  mature  evergreen  forest  can  reach  a 
height  of  at  least  120  feet.  The  thickness  of  the  closed  canopy  is  about 
40  to  50  feet,  the  same  as  the  deciduous  forest.  Thus,  the  bottom  of  the  canopy  can 
be  as  high  as  70  to  80  feet  above  the  ground.  Typical  tree  models  are  shown  in  Figure 
1(a)  and  Figure  2.  A  conifer  can  keep  its  leaves  for  three  to  eight  years.  Both  new 
and  old  leaves  are  therefore  staying  simultaneously  on  the  branches.  As  a  consequence, 
the  evergreen  tese  ho;  a  t.-uch  greura;  leaf  mass  then  o  deciduous  ree  of  the  seme  size. 
Since  the  evergreen  forest  has  a  much  greater  leaf-area,  and  the  leaves  are  functioning 
most  of  the  year  for  the  assimulation  of  bio-mass,  it  has  a  greater  growth  rate  than  a 
deciduous  forest.  For  this  reason,  management  of  evergreen  forest  crops  fer  timber  has 
long  been  in  practice.  Such  managed  forests  are  widespread  geographically. 

In  a  managed  forest,  the  tree  stands  can  either  be  even-aged,  or  graded  in  age  and 
size  by  means  of  rotation  and  selection.  Such  forests  can  achieve  their  optimum 
density  of  foliage  for  maximum  forest  yield.  It  is  important  to  note  that  in  a  young 
even-aged  conifer  stand  the  trees  have  almost  as  much  foliage  as  a  mature  conifer 
stand.  In  comparison  with  a  mature  forest,  a  young  tree  stand  has  a  higher  number 
density,  a  larger  crown  height  to  width  ratio,  and  a  larger  number  of  leaves  per 
branch.  In  fact,  leaves  and  small  limbs  of  a  young  conifer  stand  may  account  for 
50%  to  60%  of  the  total  plant  mass  above  the  ground. 

Other  than  the  major  types  of  forest  as  discussed  above,  there  are  a  iso  shrubs  and 
orchards  which  may  cover  large  areas  of  land.  In  general,  the  temperate  forests, 
including  small  trees  and  shrubs,  have  a  homogeneous  and  relatively  simple  structure. 
The  main  canopy  is  commonly  uniform  in  thickness,  and  has  clearly  defined  top  and 
bottom.  Undergrowth  density  is  sparse,  and  its  height  is  moderately  below  the  bottom 
of  the  main  canopy.  Density  of  the  canopy  depends  on  the  species  and  seasons.  From 
a  well  stocked  evergreen  forest  to  a  deciduous  forest  in  winter,  the  foliage  density 
varies  greatly.  Correspondingly,  their  effects  on  ground  attenuation  of  sound  can 
also  be  significantly  different. 

The  Tropical  Rain  Forest 

In  a  tropical  rain  forest,  the  plant  community  is  extremely  complicated.  Some  times, 

Jr  seems  to  be  such  a  chaos  that  nature  appears  to  be  anxious  to  fill  every  available 
space  with  stems  and  leaves.  Nevertheless,  order  prevails  in  the  rain  forests  upon 
close  studies  of  the  arrangement  of  the  forest  canopy  and  the  ecology  of  the  plant 
community.  Unlike  the  temperate  forest,  one  seldom  finds  less  than  forty 
species  of  plant  life  in  any  particular  tropical  rain  forest  locality.  Most  of  the 
plants  are  woody,  and  have  thick  leathery  leaves.  This  may  very  well  be  the  result 
of  physical  adaptation  to  the  humidity,  temperature,  and  the  strength  required  to 
withstand  heavy  downpours  in  the  tropica!  zone. 


Most  of  the  tropica!  rain  forests  are  evergreen  in  appearance.  However,  there  are 
also  deciduous  tree  species  in  a  tropical  rain  forest.  A  single  tree,  or  a  group  of 
trees  m  jy  suddenly  change  color  in  their  foliage,  and  lose  all  the  leaves  for  a  short 
period  of  time,  say,  from  a  few  days  to  a  few  weeks.  The  leaf  renewal  cycle 
is  not  necessarily  annual.  It  can  lie  between  eight  months  to  two  and  a  half 
years,  and  is  not  uniform  even  for  trees  of  the  same  species.  Also,  renewal  and 
shedding  of  leaves  are  not  strongly  related  to  the  dry  seasons  in  the  tropical  zone. 
Normally,  all  the  deciduous  species  are  large  trees.  Nearly  all  the  undergrowth 
plants  are  evergreen. 

In  a  tropical  rain  forest,  the  crowns  of  the  trees  are  arranged  in  many  stories. 

Although  no  standard  of  classification  has  been  agreed  upon  in  the  literature,  some 
general  description  is  possible.  According  to  Richards  (Reference  7),  most  of  the 
mixed  tropical  forest  stands  have  a  so-called  A,  B,  C  structure.  The  A -story  is 
comoosed  mainly  of  standard  trees.  These  trees  are  very  large,  and  their  crowns 
develop  freely.  However,  the  crowns  in  this  story  seldom  form  a  closed  canopy.  The 
B-story  is  composed  of  dominant  and  codominant  frees.  A  thick  and  nearly  closed 
canopy  is  formed.  There  are  also  the  smaller  trees,  not  necessarily  of  the  same  species 
as  the  dominant  or  codominant  frees,  which  will  form  a  lower  canopy,  and  is  called 
the  C-story.  A  typical  configuration  is  given  in  Figure  3. 

Other  them  the  mixed  forest,  there  are  also  some  tropical  rain  forests  which  are 
dominated  by  a  single  species  of  large  frees.  The  single-dominant  forests  have  a 
somewhat  different  structure.  Instead  of  the  A,  B,  C-structure  as  described  above, 
such  forests  can  be  said  to  have  an  A, C-structure. 

The  condition  of  undergrowth  inside  a  tropical  rain  forest  is  much  more  open  than  the 
cluttered  and  impenetrable  appearance  as  described  in  many  literary  writing.  There 
are  two  possible  reasons  that  such  impressions  were  conceived  by  earlier  writers  and 
travellers.  It  is  possible  that  most  of  the  early  travellers  made  their  passage  through 
the  tropical  rain  forests  on  waterways  and  trails.  Both  light  and  water  are  richly 
available  along  the  banks  of  streams  and  rivers.  Hence  the  undergrowth  can  expand 
to  an  impressive  thicknesfs  at  ail  heights  under  the  tree  canopies.  Similar  conditions 
prevail  along  trails  where  sunlight  is  readily  available.  Another  reason,  as  offered 
by  Richards,  is  that  the  exaggerated  manifestation  of  tropical  plant  forms  would 
tempt  the  early  explorers  to  write  accordingly.  In  general,  the  undergrowth  in 
tropical  rain  forests  is  much  denser  than  those  in  a  temperate  forest.  It  contains  also 
a  much  larger  variety  of  species  of  plant.  However,  passage  on  foot  through  the 
jungle  floor  is  not  too  difficult,  though  if  may  be  necessary  sometimes  to  remove  a 
low-hanging  branch.  Visibility  inside  the  tropical  rain  forest  is  generally  fair. 

It  depends  mainly  on  the  density  of  foliages  at  eye  level.  Hence,  visibility  is 
not  an  accurate  indicator  of  the  density  of  plant  masses  in  the  forest  as  a  whole. 

These  comments  can  be  applied  equally  well  to  visibility  conditions  in  a  temperate 
forest. 
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rimary  Mixed  Dipterocarp  Forest, 
strip  of  forest  200  ft  long  and  25 
•  \After  Richards,  Reference  7.) 


A  special  feature  in  tropical  rain  forests  is  the  abundance  in  growth  of  wood/  lianas. 
These  climbing  species  can  reach  great  heights  and  their  stems  can  have  diameters  of 
up  to  a  few  inches.  Their  interaction  with  the  tree  species  in  the  tropical  rain  forest 
can  become  an  important  part  of  the  overall  plant  ecology  in  certain  locations.  The 
abundance  of  dependent  and  serni-dependent  species  is  also  characteristic  of  the  rain 
forests.  The  presence  of  these  plants  often  adds  significantly  to  the  dramatic  appear¬ 
ance  of  a  tropical  rain  forest. 

In  general,  tropical  rain  forests  include  not  only  the  "wet"  evergreen  forests  as 
described  above,  but  also  some  other  forms  of  forests  under  limiting  conditions.  In 
parts  of  the  tropical  zone,  rainfall  and  mineral  resources  in  the  soil  can  fall  far  below 
their  normal  conditions.  Depending  on  the  severity  of  these  limitations,  a  typical 
tropical  rain  forest  can  degenerate  into  a  seasonal  evergreen  forest,  a  seasonal  semi¬ 
evergreen  forest,  a  deciduous  forest,  shrubs  and  small  trees,  or  simply  a  field  of  thorn 
and  cactus.  A  figure  showing  the  general  configuration  and  plant  density  in  these 
forms  of  tropical  forests  is  adapted  from  Reference  7  (Figure  4). 

3.3  The  Grain  Crop  and  the  Grass  Fields 


Grain  fields  and  pastures  are  perhaps  the  most  expansive  types  of  ground  cover. 

Inspite  of  the  large  variety  in  grain  and  grass  species,  their  structural  appearances  are 
generally  similar.  The  entire  canopy  is  commonly  composed  of  purely  a  single  species, 
perhaps  at  most  two  to  three;  it  is  homogeneous;  and  it  has  a  uniform  canopy  thickness. 
Most  of  the  field  crops  are  annual  plant  species:  woody  growth  is  uncommon  except 
special  cases  such  as  cotton  and  alfalfa.  Leaves  form  the  major  portion  of  the  plant 
mass  in  the  canopy. 

The  height  and  density  of  the  canopy  depend  mainly  on  the  species  and  the  practice  of 
cultivation.  Tall  species  include  corn,  wheat,  barley,  and  sugar  cane;  medium  height 
species  include  rice,  alfalfa,  and  timothy;  short  species  include  many  common  grasses 
such  as  clover  and  bluegrass.  Canopies  with  the  highest  densities  are  sugar  cane  and 
corn  for  the  tall  species,  and  rice  and  clover  for  the  short  species.  Actually,  the 
density  in  the  canopy  depends  greatly  on  the  stage  of  growth  and  the  moisture  content 
at  a  particular  time.  The  detailed  quantitative  description  of  the  structure  of  field 
crops  will  be  discussed  in  the  next  section. 
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THE  ACOUSTIC  PROPERTIES  OF  GROUND  COVER 

Materials  commonly  used  in  noise  reduction  or  noise  control  engineering  are  mostly 
porous  acoustical  materials.  Basically,  such  materials  are  composed  of  either  a  fiber 
matrix  where  air  may  pass  through  the  remaining  empty  space,  or  simply  c  homogeneous 
solid  with  numerous  open  and  close  ended  air  passages  embedded  in  it.  When  sound  is 
propagating  in  such  a  medium,  the  air  in  the  void  is  assumed  to  be  the  main  wave 
carrier.  At  the  passage  of  sound  wave,  the  friction  between  the  moving  air  particle 
and  the  surface  of  the  solid  matrix  will  cause  a  loss  of  acoustic  energy.  Moreover, 
the  elastic  properties  of  wave  propagation  in  small  air  passages  can  be  significantly 
different  from  the  corresponding  properties  in  free  air.  In  many  porous  acoustical 
materials,  heat  transfer  occurs  between  ths.  air  and  the  solid.  The  bulk  elasticity 
coefficient  of  air  will  take  its  value  under  isothermal  conditions,  instead  of  its  usual 
value  under  adiabatic  conditions  for  sound  propagation  in  free  air.  The  solid  matrix 
may,  or  may  not,  respond  to  the  sound  pressure  fluctuation.  For  some  porous  material, 
the  matrix  has  a  rigid  structure,  and  thus  remains  stationary  for  wave  passage  at  all 
frequencies.  For  others,  the  structure  of  the  matrix  is  relatively  compliant.  At  low 
frequencies,  the  solid  matrix  oscillates  with  the  surrounding  air.  Each  of  these  types 
cf  porous  acoustical  materia!  has  different  acoustical  properties. 


In  Reference  11 ,  Beranek  introduced  a  method  to  estimate  the  acoustical  properties  of 
porous  materials.  The  structural  properties  of  a  typical  porous  acoustical  material  can 
be  represented  by  the  following  parameters: 

p  =  density  of  acoustical  material 

m  ' 

Y  =  porosity,  the  ratio  of  volume  of  voids  in  material  to 

total  volume 

k  =  structure  factor,  an  empirical  constant  to  indicate  the 

nature  interstices  in  skeleton 

R  =  alternating  flow  resistance  for  unit  thickness  of  material 

Q  =  volume  coefficient  of  elasticity  of  acoustical  materia! 

p  =  density  of  air 

0 

K  =  volume  coefficient  of  elasticity  of  air  for  isothermal 

conditions. 

In  Reference  11,  all  units  are  defined  in  mks  units.  The  unit  for  alternating  flow 
resistance  is  defined  as  mks  rayls.  Two  types  of  acoustical  materials  are  considered  by 
Beranek:  materials  with  a  rigid  matrix;  and,  materials  with  a  soft  matrix  where 
K  >  20Q.  Important  acoustical  parameters  are  defined: 
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(b  +  j  b  )  -  wave  propagation  constant 
R  I 

characteristic  impedance  of  acoustical  material 
speed  of  sound  in  the  material 

proportional  to  the  real  part  of  the  wave  propagation  constant; 
it  represents  the  attenuation  coefficient  of  sound  in  the  material 
per  unit  distance 


u  =  27rf  = 


frequency;  radiar/sec 


n  =  refraction  coefficient 

j3  =  specific  admittance  ratio  with  respect  to  air 

c  =  speed  of  sound  in  air 

0 

Z  =  characteristic  impedance  of  air 

0 

The  wave  propagation  constant  is  defined  such  that  the  sound  pressure  fluctuation  can 
be  represented  as  a  function  of  distance: 

p(x)  =  exp  (-bx) 

Generally,  the  sound  wave  propagates  in  the  porous  material  at  a  speed  lower  than  the 
speed  of  sound  in  air,  and  it  attenuates  as  an  exponential  function  of  disrance  or  time. 
For  materials  with  a  rigid  matrix,  the  following  equations  are  given: 


p  kY 

0  _ 

K 


p  k  u 
0 


Z  - 
m 


.  bK 


For  materials  with  a  soft  matrix,  where  K  >  20Q  ,  the  equation  for  the  wave 
propagation  constant  is: 


b  =  ju  — 


<P,>  -j 


<  R  > 
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where 


<  R  > 

i 


R,  { 

1  - 

p.(’-Y>/pj 

{w  M)/pm : 

to  m  . 

h 

[  1  +  (R/p^u)2  1 1  +  po  (k-1 )/  pm  |'2} 

02) 


and 


<  P  >  =  P  k 
1  0 


1  t(R,/pm“)2  (Y  4  pm/p»  k  )  {'  *  po  M)/pm( 
'  +  (R,/pm“)2{'+Po  <M)/pm}'J 


(13) 


The  quantities  <  y  and  <  >  are  called  the  effective  density  and  rhe  effective 

flow  resistance,  respectively.  Other  acoustical  parameters  can  be  derived  from  the 

values  of  b  and  Z  : 

m 


c  =  t/b  ;  n  =  c  / c  ;  3  =  Z  /Z 
m  I  o  m  )  o  it* 


(14) 


Beranek  had  used  these  formulas  to  predict  the  acoustical  parameters  of  a  number  of 
commercial  acoustical  porous  material.  The  results  agreed  reasonably  well  with 
experimental  measurements. 

The  structure  of  a  forest  canopy  is  typical  of  a  porous  acoustical  material.  The  solid 
matrix  is  represented  by  stems,  branches,  and  the  foliage  of  the  plants.  The  total 
area  of  exposed  surface  area  of  leaves  and  stems  is  appreciable,  and  a  significant 
value  of  alternating  flow  resistance  :s  expected. 

The  acoustical  properties  of  ground  cover  can  be  estimated  by  using  the  Beranek 
equations.  It  is  necessary,  however,  first  to  establish  the  structural  parameters  such 
as  porosity,  effective  matrix  density,  and  flow  resistance.  The  evaluations  of  these 
parameters  are  discussed  separately  as  given  below. 

4.1  Porosity 

In  a  natural  forest  stand,  the  growth  of  the  plants  in  the  community  is  governed  by  the 
availability  of  light,  mineral,  water,  and  by  their  interdependence  and  ihe  process  of 
selection.  The  maximum  amount  of  plant  masses  which  can  survive  over  a  unit  area  of 
land  is  limited  by  these  factors.  In  a  well  stocked  forest  stand,  the  average  distance 
between  two  neighboring  trees  is  approximately  12,5  times  the  average  tree  diameter 
at  breast  height  (d.b.h.,  breast  height  is  defined  as  4.5  feet  or  1 .3  meters  above  the 
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ground).  Most  of  the  solid  volume  under  the  canopy  is  occupied  by  the  trunks  and 
branches.  The  rest  of  the  solid  material  will  be  leaves,  twigs,  lianas,  and  under¬ 
growth.  The  ratio  between  branchwood  volume  and  stemwood  volume  is  about  one  to 
four,  depending  on  the  condition  of  the  individual  foresis.  Hence,  if  the  tree  trunk 
is  assumed  to  be  conical  in  shape,  it  is  easy  to  calculate  that  approximately  0.3 
percent  of  the  available  space  under  the  canopy  is  occupied  by  solid  materials. 
According  to  accurate  measurements  of  timber  volume  in  various  forest  stands,  the 
solid  to  total  volume  ratio  can  be  as  high  as  0.45  to  0.52  percent. 

On  the  other  hand,  some  mixed  deciduous  forest  stands,  such  as  the  oak-chestnut 
forest  in  the  southeast  United  States,  the  average  tree  distance  to  diameter  ratio  can 
be  as  large  as  20.  In  such  cases,  the  solid  volume  to  total  volume  ratio  can  be  as  low 
as  0.1%.  In  Table  1,  some  typical  statistics  for  forest  stands  of  various  species  is 
presented.  In  a  forest,  the  total  quantity  of  solid  material  represents  the  accumulated 
yield  from  biological  activities  over  a  period  of  many  years.  For  the  canopies  of 
grain  crops  and  pastures,  the  conditions  are  different.  The  total  fresh  weight  of  the 
plants  is  commonly  between  5  to  15  tor/acre,  depending  on  the  canopy  height  and 
plant  density.  Hence,  the  solid  volume  to  total  volume  ratio  is  limited  to  below 
approximately  0.3%  .  In  many  cases,  this  ratio  can  be  as  low  as  0.05  to  0.1  %  . 

Hence,  the  porosity  of  a  natural  ground  cover  can  have  a  value  somewhere  between 
0.995  and  0.999. 

4.2  Effective  Density 

From  Equations  (9)  to  (13),  it  can  be  seen  that  the  density  of  the  matrix  materia!  is 
meaningful  only  if  such  masses  can  respond  to  the  sound  pressure  fluctuations  in  the 
air.  It  has  been  established  in  a  study  by  Embleton  (Reference  12)  that  even  very 
small  branches  respond  only  slightly  to  sound.  Hence,  it  is  reasonable  to  assume  that 
only  the  foliage  in  the  canopy  will  respond  to  sound  pressure  fluctuation  at  very  low 
frequencies,  following  the  definition  of  Q  in  the  Beranek  model  of  porous  acoustical 
material.  In  the  literature,  leaf  mass  has  been  measured  in  various  experimental  studies 
in  forestry  (References  5  through  10). 

It  is  important  to  note  that  the  strategy  of  leaf  growth  on  any  plant  is  governed  perhaps 
entirely  by  the  availability  of  sunlight.  The  photosynthesis  process  is  fully  effective 
only  if  the  available  light  intensity  is  above  20%  of  the  average  sunlight  intensity. 

In  permanent  shades  where  the  light  intensity  falls  much  below  this  value,  plant 
growth  will  be  impossible  except  for  a  few  species  of  plant  with  extremely  high  shade 
tolerance.  It  is  explained  by  Horn  (Reference  9)  that  there  can  be  two  strategies  for 
leaf  arrangement  for  plan*  growth  over  a  fixed  unit  area  of  land:  the  multilayer 
strategy  and  the  monolayer  strategy.  For  an  object  in  the  sunlight,  a  shadow  will  be 
casted  behind  it.  However,  the  "hard"  shadow  will  disappear  after  30  or  40  diameters 
away,  because  the  sun  is  a  light  source  of  extremely  large  dimension.  If  a  tree  has 
adopted  the  multilayer  strategy,  then  the  leaves  on  the  top  layer  will  have  an  area 
density  much  less  than  the  ground  surface  area,  say,  50% .  About  40  to  50  average 
leave  diameters  away,  a  second  layer  of  leaves  will  be  grown.  Since  the  sunlight 
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TABLE  1 


DENSITY  OF  VARIOUS  TREE  SPECIES  FOR  AVERAGE  SITES 
AT  AGES  OF  CULMINATION 


Species 


Average 
b.h.,  inch 

Trees  Per  Acre 

Basal  Area  Per 
Acre,  sq  ft 

Average  Separa 
Distance  to 
d.b.h.  *  Rati 

14.3 

360 

437 

8.9 

10.4 

628 

372 

9.6 

3.3 

3700 

220 

12.5 

10.0 

585 

317 

10.4 

5.0 

1600 

195 

13.3 

6.1 

1130 

216 

12.6 

1.8 

7400 

140 

15.6 

5.8 

1190 

221 

12.5 

4.1 

1800 

162 

14.6 

4.5 

1480 

158 

14.7 

4.9 

1090 

148 

15.2 

7.0 

540 

144 

15.4 

3.8 

1490 

118 

17.1 

5.4 

800 

122 

16.8 

3.1 

1900 

100 

18.5 

3.4 

1680 

108 

17.8 

3.8 

1150 

93 

19.2 

8.5 

320 

126 

16.5 

2.5 

2240 

68 

22.4 

4.0 

965 

84 

20.2 

6.7 

390 

90 

19.5 

Red  Fir 
Redwood 
Red  Gum 
White  Fir 
White  Pine  (Wis.) 
Sitka  Spruce 
S.  White  Cedar 
W. White  Pine 
Red  Spruce 
Shortleaf  Pine 
Slash  Pine 
Loblolly  Pine 
Lodgepole  Pine 
Douglas  Fir  (NW) 
Ponderosa  Pine 
Jack  Pine 
Long  leaf  Pine 
E.  Cottonwood 
Virginia  Pine 
Oak 

(Central  States) 
N.  Hardwood 

(Lake  States) 


*  d.b.h.  =  Diameter  at  Breast  Height;  Breast  Height  is  defined  as 
4.5  ft  or  1.3  meters  above  the  ground. 
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coming  through  the  gaps  between  the  leaves  in  tl.e  first  layer  will  cast  nc  hard  shadow 
at  the  location  of  the  second  layer,  the  sunlight  will  appear  as  a  diffused  light  source 
with  reduced  intensity.  By  repeating  the  same  tactics,  more  layers  of  leaves  can  be 
grown  profitably  until  the  diffused  sunlight  intensity  falls  below  the  20%  limit.  In 
the  monolayer  strategy,  a  plant  will  grow  as  much  non-overlapping  leaves  as  possible 
in  one  layer.  All  the  available  sunlight  will  be  intercepted  by  this  layer,  and  a  hard 
shadow  will  be  casted  on  the  ground.  Leaf  growth  will  not  be  possible  below  the  mono- 
layer  of  leaves. 

Mathematically,  the  multilayer  strategy  can  utilize  the  available  sunlight  more 
efficiently.  Because  the  effective  interception  of  light  in  consecutive  layers  forms  a 
geometrical  progression  with  a  ratio  of  at  least  1/2,  the  maximum  leaf  area  per  unit 
land  area  is  limited.  In  nature,  it  is  found  that  the  multilayer  strategy  is  generally 
followed  by  large  trees  and  plants  of  higher  ranking  species.  However,  for  the  under¬ 
growth  in  a  forest,  seedling  trees,  and  other  small  plants,  the  monolayer  strategy  will 
gather  light  most  efficiently  for  an  individual  plant  because  very  low  average  light 
intensity  prevails  in  such  an  environment.  In  Table  2,  the  leaf  area  and  leaf  mass 
for  various  important  species  of  trees  have  been  summarized.  Most  of  the  given  values 
are  obtained  from  the  forestry  literature  by  direct  or  indirect  estimates.  Their  accuracy 
is  only  nominal.  In  Table  3,  typical  values  of  leaf  area  to  ground  area  ratio  for  grains 
and  grasses  are  tabulated.  The  average  weight  per  unit  area  of  grains  and  grasses  has 
not  been  tabulated  because  the  value  varies  greatly  with  plant  density,  moisture,  and 
seasonal  conditions.  The  overall  value  for  the  common  field  crops  is  estimated  as 
approximately  5  to  15  ton/acre  of  fresh  weight. 

4.3  Alternating  Flow  Resistance 


Direct  measured  values  of  alternating  flow  resistance  of  plants  and  foliages  are  not 

available  in  the  literature.  An  attempt  is  made  here  to  obtain  indirect  estimates  of 

such  values.  As  air  flows  through  a  porous  acoustical  material,  the  boundary  layer 

buildup  at  the  exposed  surfaces  of  the  matrix  material  provides  the  mechanism  for 

momentum  dissipation.  In  an  ordinary  commercial  acoustical  material,  the  exposed 

surface  area  per  unit  volume  is  in  the  order  of  100  to  10,000  ft-1  .  In  Beranek 

(Reference  11),  the  alternating  flow  resistance  for  various  types  of  acoustical 

material  at  various  densities  has  been  measured.  It  appears  from  these  results  of 

measurement  that  the  functional  dependence  of  the  alternating  flow  resistance  on 

density  is  practically  the  same  for  all  materials.  On  a  double-logarithmic  scale, 

all  the  R  versus  p  curves  are  parallel  to  each  other.  Their  difference  in  point  of 
1  m 

origin  is  strongly  related  to  the  characteristic  value  of  exposed  area  per  unit  mass  for 
each  given  material.  Some  of  the  results  given  in  Reference  11  are  extrapolated  for 
lower  densities.  The  curves  are  shown  in  Figure  5.  By  using  the  nominal  values  of 
fiber  diameter,  the  exposed  area  per  unit  solid  mass  has  been  estimated  for  the  Aerocor 
Fiberglas  materials  at  a  density  of  0.5  Ib/cu  ft,  the  minimum  surface  area  per  unit 
volume  is  about  112  ft-1  ,  Thus,  the  alternating  flow  resistance  can  be  replotted  as  a 
function  of  exposed  surface  area  per  unit  volume  (see  Figure  6).  In  Reference  11 ,  data 
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LEAF  WEIGHT  AND  LEAF  AREA  FOR  VARIOUS  FOREST  TREE  SPECIES 
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For  broadleaf  trees,  the  area  is  measured  only  on  one  side  of  the  leaf 


TABLE  3 


I 

I 

I 


LEAF  AREA  TO  GROUND  COVER  AREA  RATIO 
FOR  GRAINS  AND  GRASSES 


Species 

Leaf  Area  Ratio  * 

Proctor  Barley 

8.1 

Winter  Wheat 

6.2 

Spring  Wheat 

4.1 

Herta  Barley 

3.7 

Rice 

10.2 

Sugar  Cane 

11.8 

Corn 

8.5 

Clover  Grass 

5.4 

Reed  Canary  Grass 

3.7 

*  The  area  is  measured  on  only  one  side  of  the  leaf. 


i 

1 

1 
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Effective  Density  of  Acoustic  Material  -  kg/m3 

Figure  5.  Air  Flow  Resistance  as  c  Function  of  Acoustic  Material  Density. 

Materials  shown  in  this  figure  are:  1  -  PF  105  XAA  Fiberglas  (1950); 

2  -Aerocor  Fiberglas  (1954);  3  -  Johns  Manville  Spintex  (1954);  and 

4  -  Canopy  of  Forest,  Grain  Crop,  and  Grass. 


Coarse  Wire  Mesh  Screens 
(Reference  1 1 ) 


has  also  been  given  for  the  flow  resistance  of  coarse  wire  screen  meshes.  If  one  assumes 
that  the  exposed  surface  area  of  the  wire  mesh  were  distributed  not  at  the  plane  of  the 
mesh,  but  throughout  a  layer  of  unit  thickness,  then  the  corresponding  flow  resistance 
value  can  be  represented  as  a  data  point  in  Figure  6.  The  value  of  flow  resistance  for 
the  coarse  wire  mesh  is  chosen  because  the  distance  between  wires  is  significantly 
large  such  that  the  boundary  layers  on  neighboring  wires  are  independent  of  each 
other.  By  drawing  a  line  through  these  data  points  and  parallel  to  the  alternating 
flow  resistance  curve  of  the  Aerocor  fiberglas,  another  estimate  of  the  flow  resistance 
can  be  established.  It  should  be  noted  that  the  flow  resistance  is  proportional  to  the 
1 .34  power  of  the  value  of  surface  area  per  unit  volume.  In  acoustical  materials 
commonly  used  for  noise  control,  both  the  fiber  size  and  the  dimension  of  the  air 
passage  are  relatively  small.  The  boundary  layers  on  neighboring  fibers  merge  rapidly. 
Thus,  the  friction  loss  owing  to  boundary  effects  would  resemble  those  of  pipe  flows. 
However,  for  surface  areas  with  very  low  densities,  the  boundary  layers  on  different 
surfaces  may  remain  independent  of  each  other  at  all  times.  In  this  case,  the  flow 
resistance  should  be  directly  proportional  to  the  surface  area  in  a  unit  volume.  This 
postulation  is  represented  by  the  dash  line  in  Figure  6.  Within  the  range  of  interest 
of  the  present  study,  the  differences  among  the  above  three  estimates  are  rather  small . 

According  to  the  values  of  leaf  area  ratio  and  canopy  thickness  for  various  plant 
species,  the  leaf  area  per  unit  volume  can  be  estimated.  The  typical  ranges  of  values 
are  given  in  Figure  7.  In  some  cases,  the  exposed  surface  area  for  stems  and  branches 
may  account  for  a  significant  portion  of  the  total  surface  area  in  the  canopy.  It 
depends  on  the  total  volume  of  wood  and  the  average  branch  and  stem  sizes.  Such 
dependences  are  given  in  Figure  8.  By  using  Figures  6  through  8,  one  can  readily 
estimate  the  alternating  flow  resistance  of  a  given  type  of  plant  canopy. 

.4  Other  Parameters  and  Considerations 

There  are  three  more  parameters  remaining  to  be  considered:  the  volume  coefficient 
of  elasticity  of  the  porous  acoustical  material,  Q;  the  volume  coefficient  of  elasticity 
of  air,  K;  and,  the  structure  factor,  k.  Since  the  porosity  is  extremely  close  to  one, 
the  value  of  k  can  be  assumed  to  be  1  .0  for  all  cases.  It  is  not  necessary  to  know  the 
exact  value  of  Q  in  the  Beranek  equations.  Since  both  the  material  density  and  the 
flow  resistance  are  provided  by  the  leaves,  the  plant  canopy  is  assumed  to  be  a  soft 
porous  acoustical  material,  i.e.,  K  >  20Q.  For  the  forest  canopies,  the  foliage 
density  is  relatively  low.  The  parameter  K  ;s  assumed  to  take  its  value  under  isen- 
tropic  conditions.  For  dense  plantation  of  grain  crops  or  grasses,  K  can  possibly 
take  its  value  somewhere  between  isothermal  limit  and  isentropic  limit.  Further 
studies  are  definitely  required. 

According  to  the  above  defined  ranges  of  structural  parameters,  the  acoustic  properties 
of  ground  cover  materials  have  oeen  estimated.  Computations  are  made  assuming  both 
conditions  where  the  matrix  structure  can  either  be  rigid  or  be  soft.  The  obtained  values 
for  the  refraction  coefficient,  n,  and  the  internal  transmission  loss  coefficient,  a  , 
are  plotted  as  a  function  of  frequency. 


The  results  for  the  soft  matrix  materials  are  given  in  Figures  9  and  iO.  According  to 
the  Beranek  model  of  acoustical  materials  with  a  soft  blanket,  the  velocity  of  sound 
in  the  low  frequenc/  range  is  greatly  reduced  by  the  effective  material  density.  For 
higher  frequencies,  the  soundwave  can  propagate  in  the  acoustical  material  at 
nearly  the  speed  of  sound  in  free  air.  For  the  attenuation  coefficient  a  ,  (he  trend 
with  frequency  is  just  the  opposite.  Since  the  matrix  rides  with  the  aii  flow  at  low 
frequencies,  very  little  friction  is  generated  through  friction  between  the  air  flow 
and  the  matrix.  Hence,  a  is  small  in  the  low  frequency  range.  As  the  frequency 
increases,  the  inertial  forces  will  keep  the  matrix  material  stationary,  and  thus  the 
relative  motion  between  the  matrix  and  air  increases.  Consequently,  a  increases 
with  frequency,  and  reaches  a  limiting  value  for  very  high  frequencies. 

Values  of  the  refraction  coefficient,  n,  for  materials  with  a  rigid  matrix  are  plotted 

in  Figure  11 .  For  a  rigid  matrix  acoustical  material,  the  acoustical  properties  do  not 

depend  on  p  .  Hence,  R  is  the  only  sianificcnt  variable.  The  dependence  of  a- 
m  i 

on  frequency  is  very  weak  for  acoustical  materials  with  a  rigid  matrix.  It  remains 
about  constant  for  all  frequencies.  The  value  of  a  for  an  acoustical  material  with 
a  rigid  matrix  is  approximately  the  same  as  the  limiting  value  of  a  at  high  frequencies 

for  a  soft  matrix  material  with  the  same  value  of  R  . 
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Figure  9.  Refraction  Coefficient  Versus  Frequency.  Each  curve  represents  a  soft  acoustic  material 
with  given  values  of  effective  density  and  air  flow  resistance. 


1—1 


Figure  11 .  Refraction  Coefficient  Versus  Frequency.  Each  curve  represents  a  rigid  frame 
acoustical  material  with  a  given  air  flow  resistance. 


5.0  PREDICTION  OF  ABSORPTION  SPECTRA  FOR  VARIOUS  TYPES  OF  NATURAL 
GROUND  COVER 

In  the  previous  section,  the  acoustic  properties  of  canopies  of  trees,  grain  crops,  and 
grass  fields  have  been  estimated.  According  to  these  values  of  refraction  coefficient, 
acoustic  admittance  ratio,  and  the  internal  transmission  loss  coefficient,  the  ground 
absorption  spectra  can  be  calculated. 

The  structural  properties  of  natural  ground  covers,  and  consequently  their  acoustic 
parameters,  fall  into  three  recognizable  ranges: 

•  Trees  and  shrubs 

•  Grain  crops  and  tall  grasses 

•  Short,  dense  grasses 

For  trees  and  shrubs,  the  fresh  leave  weight  per  acre  of  ground  area  is  between  4.0  to 
19.0  tons  per  acre  (Table  2).  The  leaf  area  to  ground  area  ratio  is  generally  in  the 
range  of  5  to  25.  For  tall  matured  forests,  the  effective  thickness  of  the  closed  canopy 
is  approximately  50  ft.  Therefore,  the  effective  density  for  acoustic  purposes  is  in  the 
order  of  0.075  to  0.200  kg/m3;  and  the  leaf  area  per  unit  volume  is  in  the  order  of 
0.1  to  0.5  ft-1  .  For  young  tree  stands,  tree  species  which  are  short  by  nature,  and 
shrubs,  the  effective  thickness  of  the  canopy  can  be  10  to  20  ft.  At  the  same  time, 
the  leaf  mass  and  the  leaf  area  ratio  remain  practically  the  same  as  the  tall  forests. 

Hence  the  effective  leaf  density  is  in  the  order  of  0.20  to  0.6  kg/m3;  and  the  leaf 
area  per  unit  volume  is  in  the  order  of  0.30  to  2.5  ft*1  .  In  the  summer,  the  leaf  area 
dominate  the  total  surface  area  in  the  canopy  of  a  deciduous  forest.  In  the  fall  or 
winter,  all  the  leaves  will  be  fallen.  The  only  exposed  surface  area  will  be  those  of 
the  stems  and  branches.  For  a  forest  canopy  with  a  large  portion  of  small  branches,  the 
total  surface  area  can  be  as  high  as  0.4  ft*  . 

For  grain  crops  and  tall  grass  pastures,  the  structural  parameters  fall  into  an  entirely 
different  range.  The  canopy  height  is  commonly  between  4  to  8  ft.  The  fresh  weight 
of  the  leafy  portion  of  the  plant  is  between  5  to  8  tons  per  acre,  depending  on  the 
species  and  the  plant  number  density.  The  leaf  density  is  approximately  0.4M.O  kg/m3. 
Since  all  the  leaf  shapes  in  grains  and  grasses  are  flat,  the  effective  surface  area  for 
acoustic  resistance  should  be  counted  on  both  sides  of  the  leaves.  Hence  the  effective 
leaf  area  ratio  is  twice  that  which  is  quoted  in  the  agricultural  literature  where  only 
one  side  of  the  leaf  is  measured  for  photosynthesis  purposes.  The  estimated  leaf  area 
per  unit  volume  is,  therefore  in  the  order  of  1 .0  to  4.0  ft*1  . 

The  highest  foliage  mass  density  and  leaf  area  per  unit  volume  can  be  found  in  grass 
fields  with  an  overall  thickness  of  less  than  24  inches.  The  effective  density  is  in  the 
order  of  1.0  to  2.0  kg/m  .  The  two-sided  leaf  area  ratio  is  approximately  10.  Hence, 
the  leaf  area  per  unit  volume  has  a  value  between  5.0  and  10.0  ft-1  . 
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The  ranges  of  parameters  for  the  above  three  classes  of  vegetation  are  summarized  in 
Table  4. 

Computer  programs  have  been  written  to  perform  both  the  calculation  of  acoustical 
properties  of  the  ground  cover  and  the  computation  of  sound  absorption  spectrum  using 
the  analysis  given  in  Section  2.0.  In  the  computation  of  the  sound  absorption  spectrum, 
three  acoustical  parameters  are  required:  the  refraction  coefficient,  n  ;  the  specific 
admittance  ratio  of  the  ground  cover  relative  to  air,  3  ;  and  the  normal  admittance 

ratio  of  the  ground  si rface,  3  .  The  parameters  n  and  3  are  generated  by  the 

2  1 

computer  program  according  to  the  given  structural  properties  of  the  ground  cover. 

Two  values  of  3  are  assumed  arbitrarily.  For  a  forest  stand,  3  is  assumed  to  be 
2  2 
0,5  if  the  undergrowth  is  relatively  dense;  it  is  assumed  to  be  0.2  if  the  undergrowth 
is  sparse  and  the  floor  of  the  forest  is  covered  with  only  a  shallow  layer  of  fallen 
leaves  and  branches.  For  grain  and  grass  fields,  the  ground  surface  is  expected  to  be 
a  loose  top  soil  in  most  cases.  Such  a  surface  has  a  relatively  high  acoustic 

impedance.  Hence,  the  value  of  3  is  taken  to  be  0.2. 

2 

For  forest  canopies,  the  ground  absorption  spectra  for  three  sets  of  representative 
acoustic  parameters  have  been  computed.  The  effective  canopy  thickness,  instead 
of  the  overall  height  of  the  forest  itself,  is  chosen  to  represent  the  layer  thickness 
of  the  acoustic  material  as  defined  in  the  theory.  In  these  computations,  the  receiver 
is  located  either  above  the  canopy,  or  immediately  on  top  of  the  canopy.  Computation 
has  not  been  made  for  receiver  locations  in  the  ccnopy,  although  in  practical  situations 
the  observer  may  very  well  be  standing  on  the  floor  of  a  forest  and  listening  to  the 
noise  produced  by  a  low-flying  aircraft.  Under  such  practical  conditions,  one  can 
safely  assume  that  the  minimum  excess  attenuation  is  represented  by  the  values  predicted 
for  an  observer  located  immediately  on  top  of  the  forest  canopy. 

The  computed  results  are  given  in  Figures  12  through  14,  and  Tables  5  through  19.  For 
tree  canopies  with  very  low  density  and  air  flow  resistance,  the  acoustical  effect  of  the 
canopy  on  wave  propagation  is  obviously  very  low.  In  Figure  12,  the  absorption 
spec'.um  for  small  canopy  thicknesses  exhibits  clearly  the  characteristics  of  interference 
effect  for  wave  reflection  next  to  a  solid  wall.  It  appears  also  that  the  tree  formations 
offer  very  lirtle  attenuation  to  sound  at  higher  frequencies,  owing  to  the  ground  effect. 
The  attenuation  to  higher  frequencies  is  significant  only  if  the  sound  propagates  through 
the  leafy  canopy  itself.  Significantly  large  attenuation  is  indicated  by  the  transmission 
loss  coefficient,  a  ,  as  given  in  Figure  10.  In  the  winter  time,  perhaps  very  little 
ground  attenuation  effect  can  be  expected  from  a  deciduous  forest  stand. 

The  absorption  spectrum  has  also  been  computed  for  several  densities  and  thicknesses  of 
grain  and  grass  canopies.  The  results  are  given  in  Figures  15-19  and  Tables  20-29. 

The  range  distances  and  sound  sources  height  have  been  limited  to  values  at  which  field 
measurements  are  commcnly  conducted.  These  tables  and  graphs  are  considered  to  be 
reasonably  accurate.  They  can  be  employed  for  preliminary  estimates  of  ground 
absorption  of  sound  in  actual  field  conditions. 
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TABLE  4 


MECHANICAL  PROPERTIES  OF  THREE  CLASSES 
OF  PLANT  CANOPIES 


Type  of 

Foliage  Density 

Surface  Area  Per  Unit  Volume 

Air  Flow  Resistance 

Ground  Cover 

kg/m 

ft"' 

mks  rayl/m 

Trees  &  Shrubs 

Tall  Canopy 

0.075  -  0.200 

0.1  -0.50 

1. 0-5.0 

Short  Canopy 

0.20  -0.60 

0.3  -2.50 

3.0-28 

Stems  &  Branches 

0.0 

0.05  -0.40 

0.0 -4.0 

Grain  Crops 
&  Grasses 

0.40  -1.00 

1.00-4.0 

8.0-50 

Dense  Short  Grass 

1.0  -2.0 

5.0  -10.0 

50-120 
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If  the  source  he'nht  were  different  from  those  as  given  in  the  computed  conditions, 
the  excess  attenuation  can  be  adjusted  by  using  the  formula 


E  =  E0  -20log10  <H/H() 

where  E  is  the  excess  attenuation  corresponding  to  a  height  of  H  ,  and  E  denotes 
0  0 

the  required  excess  attenuation  corresponding  to  a  height  of  H.  However,  the  above 

equation  is  accurate  only  if  both  H  and  H  are  at  least  one  layer  thick  above  the 

top  of  the  ground  cover.  In  the  tabulated  results  of  the  absorption  spectrum,  the 
excess  attenuation  is  given  at  various  range  distances.  It  can  be  noted  from  these 
results  that  the  excess  attenuation  increases  at  6  dB  per  doubling  the  range  distance 
in  the  far-field,  and  increases  at  a  somewhat  smaller  rate  in  small  range  distances 
from  the  sound  source. 


6.0  CONCLUSIONS  AND  DISCUSSIONS 


There  are  certain  limitations  to  the  present  prediction  scheme.  First  of  all,  the 
Beranek  model  of  porous  acoustical  material  is  a  reliable  method  for  the  estimates  of 
acoustical  properties  of  the  ground  cover,  yet  it  is  not  the  only  available  approach, 
for  example,  in  the  high  frequency  range,  the  transmission  loss  of  sound  through  a 
plant  canopy  may  depend  on  the  diffraction  effect  of  leaf  surfaces.  A  recent  study  by 
Ay  lor  (Reference  13)  has  shown  remarkable  agreement  between  theory  and  experiment 
for  sound  transmission  losses  through  high  density  brocdleaf  canopies.  In  the  approach 
taken  by  Aylor,  diffraction  and  penetration  of  sound  through  the  leaf  mass  are  consider¬ 
ed  to  be  chief  mechanisms.  The  improvement  of  the  definition  of  acoustical  properties  for 
ground  cover  material  can  be  achieved  by  both  theoretical  and  experimental  approaches. 
The  experiments  may  include  the  measurement  of  physical  properties  of  leafy  canopies, 
from  which  the  acoustical  properties  can  be  estimated,  or  the  direct  measurement  of 
the  acoustical  properties  of  the  plant  canopy.  The  improved  knowledge  of  ground  cover 
acoustical  properties  can  definitely  be  used  advantageously  in  the  present  layered  media 
approach  of  dealing  with  the  ground  absorption  of  sound. 

From  the  computations  of  Section  4.0,  it  is  obvious  that  ail  the  ground  cover  materials 
have  very  small  acoustical  "density"  as  compared  to  ordinary  architectural  acoustical 
materials.  However,  the  lack  of  acoustical  "density"  is  somewhat  compensated  by  the 
large  spatial  dimensions  often  encountered  in  field  conditions.  From  the  computed 
results,  it  becomes  apparent  that  the  ground  absorption  effect  is  significant  only  if  the 
range  distance  is  larger  than  approximately  25  times  that  of  the  height  of  the  sound 
source.  At  range  distance  of  less  than  10  times  that  of  the  height  of  the  sound  source, 
net  only  the  absorption  is  small  in  value,  but  also  the  spectrum  is  irregular  and  exhibits 
no  apparent  trend.  It  is  not  clear  whether  the  irregularity  is  caused  partly  by  nature  of 
the  ground  absorption  effect  or  that  it  is  caused  entirely  by  the  inaccuracy  of  the 
mathematical  model  in  the  near  field.  Fortunately,  such  near  field  conditions  are 
seldom  of  practical  interest  for  noise  control  studies. 

Inspite  of  the  limitations  as  mentioned  above,  the  computed  results  seem  to  agree  well 
with  experimental  evidences.  It  was  pointed  out  in  the  earlier  publication  under  this 
program  (Reference  1 )  that  the  predicted  ground  absorption  effects  agreed  closely  with 
the  result  of  laboratory-scaled  experiments.  In  Section  5.0,  there  are  several  sets  of 
computations  where  the  input  parameters  are  typical  of  dense  grasses  and  low  bushes. 

The  height  of  the  source  and  the  receiver  is  assumed  to  be  between  5  to  1 0  ft;  the 
thickness  of  the  ground  cover  is  from  6  in.  to  6  ft;  and  the  range  distances  are  250  and 
500  ft.  The  spectrum  shapes  in  the  predicted  cases  are  strikingly  similar  to  those 
measured  in  the  field.  The  peak  is  located  between  100  and  400  Hz,  and  the  excess 
attenuation  is  close  to  15  dB  for  all  cases.  Similar  values  of  ground  absorption  under 
similar  conditions  have  been  observed  before  (References  14  and  15).  It  is  important- 
to  note  also  that  the  transmission  loss  coefficient  for  wave  propagation  through  the 
canopy  itself  has  been  computed  for  a  range  of  values  in  Section  4.0.  For  conditions 
typical  in  a  tree  stand,  =  0,225  kg/m3,  and  R  =  20  rayl/m.  The  computed 
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transmission  loss  coefficient  is  about  5.5  dB/100  ft.  This  value  is  very  close  to  those 
observed  previously  by  various  authors  (References  12,  16  and  17). 

Several  important  trends  have  also  been  observed  in  this  study.  The  peak  of  the 
absorption  spectrum  is  governed  by  several  principal  factors:  the  height  of  the  sound 
source*;  the  thickness  of  the  ground  cover  layer;  and,  the  acoustical  properties  of 
the  ground  cover  and  the  ground  surface.  In  previous  computations  by  Ingard,  the 
absorption  peak  is  normalized  with  respect  to  the  sound  source  height.  The  product  of 
the  wavenumber  and  the  source  height,  kh,  has  a  value  of  approximately  1 1 .0.  In 
the  present  study,  the  peak  absorption  frequency  can  be  normalized  against  the  layer 
thickness.  It  is  found  in  the  computations  that  the  peak  of  the  ground  absorption 
frequency  is  dominated  by  the  layer  thickness,  provided  that  the  ground  cover  is 
sufficiently  dense.  The  value  of  kh  is  between  3.0  to  10.0.  For  layers  with  very 
small  density,  such  as  the  density  of  a  deciduous  forest,  the  influence  of  the  height 
of  sound  source  on  the  peak  absorption  frequency  will  remain  significant. 

The  predictions  of  ground  absorption  of  sound  over  ground  cover  with  very  large 
thicknesses,  such  as  forest  stands,  remain  to  be  compared  with  experiment  measurement. 
However,  it  is  reasonable  to  assume  here  that  the  predictions  as  given  in  Section  5.0 
should  serve  well  as  an  indicator  of  the  expected  ground  absorption  levels  in  actual 
field  conditions. 


According  to  the  principal  of  reciprocity,  the  value  of  excess  attenuation  remains 
invariant  if  the  positions  of  the  sound  source  and  the  receiver  are  exchanged.  Hence, 
it  is  not  necessary  to  mention  both  the  source  height  and  the  receiver  height  in  the 
discussions. 
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Figure  13.  Ground  Absorption  Spectrum  for  Tree  Stands  at  Various 
Canopy  Heights.  Heavy  undergrowth  is  assumed; 

32  =  0.5;  p  =  0.225;  R  =  30.  The  receiver  is 
immediately  on  top  of  the  canopy. 
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Figure  14,  Ground  Absorption  Spectrum  for  Tree  Stands  at  Various 
Canopy  Heights.  Sparse  undergrowth  is  assumed; 
i32  =  0.2;  p  -  0.225;  R  =20.  The  receiver  is 
immediately  on  top  of  the  canopy. 
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Figure  15.  Ground  Absorption  Spectrum  for  Grains  and  Tall  Grasses. 
Layer  Thickness  =  6  ft 
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Source  Height  =  1 0* 
Receiver  Height  =  10' 
Layer  Thickness  =  4' 
Range  Distance  =  500' 
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Figure  16.  Ground  Absorption  Spectr 
Layer  Thickness  =  4  ft 
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Figure  17.  Ground  Absorption  Spectrum  for  Dense,  Short  Grasses. 
Layer  Thickness  =  2  ft 
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Source  Height  =  5' 
Receiver  Height  =  5' 
Layer  Thickness  =  1 1 
Range  Distance  =  250' 
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Figure  18.  Ground  Absorption  Spectrum  for  Dense,  Short  Grasses. 
Layer  Thickness  =  1  ft 
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TABLE  6 

GROUND  ABSORPTION  SPECTRUM  FOR  TREE  STANDS  WITH  HEAVY  UNDERGROWTH 
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TABLE  7 
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12*59 

-1.31 

3.70 

9.21 

14*98 

20*67 

lb.fcb 

0*  3C 

b«34 

10.87 

16*64 

22*54 

19«95 

1.69 

6.7b 

12.29 

18.07 

23*96 

2b.  id 

2  *  8d 

7.90 

13.44 

19.22 

25*12 

31*62 

3  •  7o 

8.81 

14«34 

20.12 

26. C2 

3b.?  1 

4 * b4* 

9*53 

15*04 

20.80 

26*70 

bO  *  12 

b*  19 

10*10 

15.58 

21.33 

27.20 

3  3*10 

b  •  79 

10*60 

16.03 

21*76 

27*62 

79.49 

6.40 

11.  C9 

1 6  *  4b 

22*15 

28.01 

1  jO • C  J 

7 «  Of. 

11*60 

16*89 

82*55 

28*39 

12b. o9 

7.76 

12*17 

17.36 

82.99 

28.81 

Ibo  .‘♦O 

6*54 

12*81 

17-93 

23.50 

29*29 

19  3  *  b  j 

9.23 

13*49 

18.54 

2  4  *  Q7 

29*83 

cbl  *13 

9-59 

14*10 

19.16 

24*67 

30*43 

316.23 

9-28 

14*36 

19*60 

25*16 

30*95 

39o  *11 

8*16 

13.90 

19.51 

25*28 

31*12 

0^1*19 

6*46 

1 2  •  b4 

18.55 

24.57 

30*53 

6  jc  *  9b 

4  •  b4 

10.  b4 

16-74 

82.91 

28*98 

/3i*  •  33 

2*64 

8*28 

14.47 

20*65 

27.23 

1001*00 

0 .67 

6*03 

12*05 

18-16 

24*60 

12b£  *32 

-0*73 

3-87 

9.66 

15.69 

21*94 

1  5  8  4  •  t  J 

-2.  13 

1*87 

7.37 

13.38 

19*40 

193b » 2  o 

-3.26 

0*03 

5.20 

11-04 

16*96 

2ol 1 *2o 

-4.JG 

- 1  •  63 

3.15 

8*81 

14.69 

316c  *2/ 

•  4  «  w  7 

-3*06 

1.2b 

6.68 

12*94 

398 1*05 

-2*  8b 

-4.17 

-O.bl  ■ 

4 . 70 

10*72 

boll *86 

1 .  Go 

-4*82 

-2.09 

2.80 

8*61 

63c9  *  b4 

4.22 

-4.6b 

-3.43 

1*11 

6*57 

7943.2b 

-2.07 

-2*91 

-4.40 

-0.46 

4*79 

9999.96 

-3.  32 

0*41 

-4.79 

-1*55 

3*52 

71 


TABLE  29 

GROUND  ABSORPTION  SPECTRUM  FOR  DENSE,  SHORT  GRASSES 


S*  *  0* 137C0L+U1 

YbC=  0»?‘i0GGL  +  u0 
¥*£*  0.137C0E+G1 

n  «  0*152e0t+Co 

6  c  -  G*200COE-*-vO  *  c»GOOOOE*-OC 

5K  =  0* 100G0E  +ol 

V  *  0*996C0t+uJ 

r<  1  *  0  •  90oGG2'*-o2 

»  0  •  12(k (Jt  +i.  1  ft  =0.2 

.»  *  0*  10GoCL+^4  2 

GK  s  0 • 10GGul+L6 
HI  =  0*1 18CGE  t 


0N£- TnlKU-aCTANf.  HAND  EXCESS  ATTENUATI9N/  06, 


P  HE  J 

RANGE 

DISTANCES 

-  FEET. 

HZ 

25C 

500. 

1000* 

2000. 

40CQ. 

1  c  •  00 

-2*26 

2.77 

8.29 

14*06 

19*95 

12*69 

•  G  *  a  1 

4*25 

9.79 

15*56 

21*46 

16*85 

1  •  JO 

6  *  lo 

11*67 

17*46 

23*36 

19*95 

2*60 

7.74 

13*32 

15.13 

25*03 

26*12 

3*  33 

9*  10 

14.70 

20-51 

26*4? 

31*62 

4  *  38 

10*16 

15.76 

21.57 

27*49 

39.01 

5.  77 

10-94 

16.54 

22*35 

28-27 

5G*12 

6.3  7 

11*61 

17*09 

22*89 

28*80 

63.10 

6*04 

11-92 

17-48 

83.27 

29-17 

75,4  3 

7.24 

12*25 

17.76 

23.54 

29.43 

loG • CO 

7.62 

12*53 

18.00 

23*75 

29*64 

125.83 

7.99 

12*  8j 

18. 22 

83*95 

29.88 

150*43 

6  •  36 

1 3 » C7 

18.44 

84*13 

29*95 

193.53 

6.64 

13.32 

18*64 

24.31 

30*18 

251.13 

8.  73 

13*45 

18.76 

2-4.41 

30*30 

3lo»23 

8*43 

13*33 

18.69 

24  •  36 

30*25 

39o. 11 

7 . 60 

12*77 

18*25 

24*0? 

29*89 

631.13 

6.24 

11*66 

17*30 

23*16 

29*07 

630  *  95 

4. 5b 

1 C  *  03 

15*61 

21*72 

27*67 

/34 ,33 

2.  /? 

8. 1C 

13*92 

19*86 

26  *25 

1000*0 0 

C  *  3? 

6  •  03 

11.81 

17*75 

24.03 

123*  *  92 

- o  *  7b 

3.96 

9*62 

15.53 

21*68 

1 5  a  4 » e,  9 

-2.21 

1.92 

7.45 

13*38 

19*35 

1596.25 

-  3 »  36 

0*18 

5*34 

11*15 

17*08 

2611*33 

-  4 « 05 

-1*57 

3*32 

9  *  00 

14*87 

3162.27 

-  4 . 02 

-3*08 

1.42 

6*91 

13*18 

33a 1 . Q5 

-  C  *  uci 

-4.14 

-0.35  ' 

4.93 

10*98 

5w 1 1 ♦ 86 

1*32 

-4.77 

-1.95 

3*04 

8*90 

63Q9  *  64 

3*  94 

-4*57 

-3*30 

1.34 

6*8? 

7  9  4  3  *  2  5 

-  1  >  3  8 

-2.83 

-4.29 

-0*25 

5*06 

9939.9a 

- 3*  08 

0.3  a 

-4.71 

-1*37 

3*76 
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WYLE  LABORATORIES 
COMPUTER  PROGRAM  DESCRIPTION 


Program  Number:  T2/ 002S-I 

Author  :  D.  M.  Lister 

Date  :  January  28,  1972 

Computer  XDS  Sigma  5 

Source  Language:  Fortran  IV-H 

Monitor  System  :  RBM  -  1 

1 .0  PROGRAM  TITLE 

Computation  of  ground  attenuation  effects  on  noise  spectra. 

2.0  PURPOSE 

Given  a  mathematical  definition  of  source  height,  refraction  and  specific  admittance 
coeffients,  the  program  computes  the  spectral  sound  attenuation  for  five  (5)  horizontal 
distances  (x)  from  the  source.  The  attenuation  is  computed  ot  five  (5)  heights  for  each 
x  and  the  average  recorded.  A  third-octave  averaging  process  is  performed  in  order  to 
smooth  the  results.  The  user  has  direct  control  over  the  choice  of  spectral  and  averaged 
plots  to  be  produced . 

3.0  METHOD 

The  method  employed  is  described  in  Reference  1 . 

4.0  COMPUTER  CONFIGURATION 

The  required  hardware:  XDS  Sigma  5  computer  with  16Kcore,  card  reader,  lineprinter 
and  Calcomp  plotter. 

5.0  DATA  INPUT 

The  data  input  is  in  the  form  of  punched  cards,  the  formats  of  which  are  described  in  the 
following  table. 
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Note 

Card 

Variable 

Description 

Format 

Cord 

Columns 

1 

1 

1P2  (1) 

Plot  control  parameters  for  1/3  octave 

averaged  plots 

11 

1 

IP2  (2) 

IP2  (I)  controls  plot  for  range  f 

11 

2 

IP2  (3) 

!f  1P2  (I)  =  0  then  no  plot  for  range  1 

11 

3 

IP2  (4) 

If  1P2  (1)  =  1  then  plot  on  new  axes 

11 

4 

IP2(5) 

If  IP2  (I)  =  2  then  plot  on  existing  axes 

11 

5 

2 

Sy 

SYIN 

The  source  height 

E*1 

1-5 

YSC 

Scale  factor  for  y  “  coordinates 

6-10 

YRE 

The  receiver  height 

Hi 

11-15 

P2 

82 

The  specific  admittance  coefficient 

mm 

at  the  ground 

mm 

16-25 

HH 

Ground  cover  layer  thickness 

ESH 

26-35 

k 

SK 

Structure  factor  that  introduces  into 

equations  nature  of  intersticies  in 

skeleton 

F5.0 

36-40 

Y 

Y 

Porosity 

F5.0 

41-45 

2 

R, 

R! 

Alternating  flow  resistance 

F5.0 

46-50 

c 

RM 

Density  of  acoustical  material 

F5.0 

51-55 

3 

*  m 

Volume  coefficient  of  elasticity  of 

Q 

Q 

acoustical  material 

F5.0 

56-60 

NOTES: 

1 .  After  processing  each  set  of  data  the  program  returns  to  read  another  set  of  data  -  if  an 
END-OF-FILE  is  read  then  the  program  will  terminate. 

2.  These  variables  are  complex  numbers,  hence  the  real  part  followed  by  the  imaginary  part 
must  be  supplied. 

3.  These  variables  are  multiplied  by  1000  by  the  program,  thus  adjust  accordingly  before 
input. 

4.  See  Reference  2  for  full  description  of  this  and  following  parameters. 

5.  All  the  dimensional  units  follow  the  mks  system. 
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UNEPRINTER  OUTPUT 


6.0 

This  consist  of  an  annotated  copy  of  the  input  data,  and  a  table  of  attenuations  at 
various  range  distances  for  the  frequency  range  of  interest . 

7.0  CALCOMP  PLOTTER  OUTPUT 

This  consists  of  plots  of  the  attenuation  levels  requested  for  a  frequency  range,  on  a 
log  scale,  from  10  Hz  to  10000  Hz.  The  range  of  the  attenuation  levels  is  from  -10  db 
to  40  db  on  a  linear  scale. 

8.0  DATA  STATEMENT  OF  PROGRAM 

The  number  of  frequency  points  per  decade  (ANN)  is  preset  by  a  DATA  statement  to  be 

100. 

The  program  will  always  compute  three  decades  of  frequency  beginning  at  FRZ  which  is 
preset  to  10  Hz. 

The  spectrcl  plots  for  the  various  range  distances  may  be  obtained  by  altering  the  array 
IR  according  to  the  rules  specified  in  section  5.0  for  the  array  I  PI .  Each  member  of 
the  array  IPI  is  currently  set  to  zero. 

9.0  REFERENCES 

1.  Pao ,  S.  P.,  and  Evans,  L.  B.,  "Sound  Attenuation  Over  Simulated  Ground  Cover," 
Journal  of  the  Acoustical  Society  of  America,  Vol.  49,  pp  IQ/69-10 75,  1971, 

2.  Beranek,  L.  L.,  "Noise  Reduction,"  McGrow  Hill  Book  Company,  Inc .,  pp 257-260 
1960. 
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Figure  A-l .  Geometrical  Configuration 
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*** 


L0MPLE a  SGD,8D1 

C9MPLEX  VIA  (46)#  V2A ( 46 ) , AN1 ( 46 ) , AN2 {  46) 

COMPLEX  GZ,G1,B1,B2,SIG,SI2, Vl< VD1, V9D1, V2, vD?, VDD2,G72,G12, TEMP, R 
lZ,Rl,RZ2,Rl2,R2,K22,H,FNC.EM,EN2,Tw,TEMP2,AK,HH,ZZER 
C 

C0MM0N/L IMI TS/CZER0 
C 

CaMM0N/SW/GZ,Bl,82,  VI,  VDl,VDDl,V2,  VC2, VDD2, Gl , H, SI G, TW0PI , XZ (5  ) , Y 
1Z(5),  I  PI  <5),ANM,FKZ,  IP2(5),HH,ZZER 
C 

C0MM9N  SGD(301 ),BD1 (301 ) 

C0MM9N  FRQ(301),96P(301),X(31),Y(31,5),DB(5 ),PARM ( 7 ) , W,DUM ( 12 ) 

/— 

W 

DIMENSION  AKD ( 1 ) 

DIMENSION  hD(2),RZD(2),R1D(2),TWD(2),GZD(2),Ak(2) ,  HHD <2),B1D(2),SI 
1GD ( 2 ) , /ZERO ( 2 ) 

C 

EQUIVALENCE  (A<D,FRQ) 

EQUIVALENCE  ( hD,  h,  HR ) ,  ( HL ( 2 ) ,HI  ) ,  <R7D,RZ,R),  (RZD<2),RI  ),  (FN,FNiC),  ( 
1R1D,RR1,R1 ),  (RRI,RlD(2)  )  j  (RZ2,RR),  (RPl2,Rl2),  (  T n«  TinD,  TI*R )  ,  (TRD(2), 
2T/)I  ),  (uZC:,GZ,GZR  ),  <G7D(2),GZI  ),  { AK, A<R  ) ,  (AK(2),AKl ) 

EQUIVALENCE  ( HHD, HH,HHR ) , <HHD(2),HHI ), (BID# Bl >, <SIGD,SIG) 
lquIVAlENCE  (7ZER9,ZZER,7.ZERD),  <ZZERI,ZZERD(2>  ) 

C 

LATA  Tx0Pl/6.?«31853/ 

DATA  A, NNUFRZ^ZERe/ZZERe^ZEK  1/100*0,  10*0, 345. 31, 410. 0,0.0/ 

LATA  XZ/1215.2,6C9.6, 304.8, 152.4,76.2/ 

DATA  YZ/0. 0,0. £5,0.5,0.75, l.Q/ 

DATA  IPI  70,0,0,0,0/ 

C 

CALL  PLOT  (0.0,-12.0,25) 

CALL  PL0T  (0.0,1.5,25) 

2  F.\C»iO.  0,0*0) 

KZ2=FNC 
Kl2«F\C 
i\  R I  ■  0 » C 
N I «0«0 
bZ  I  »0  *  0 
A  <  I  *  0  •  0 
h«FNC 


C*** 

r 

READ  ( 105,3,£ND»999  )  ( IP2( I ), I =1,5) 

3  F0RMAT  (511) 

READ  (105,4)  SYIN,YSC,  YRE,B2,HH,  (PARflf  I  ),  1*1,5) 
h  F0RMAT  (16F5.0) 

FARM ( 3 ) »PARM ( 3 ) #1000*  0 
RARM(5)«PARM(5)*10G0»0 
FARM ( 6 ) »1000CO*  0 
FARM ( 7 ) »1 » 18 

rtR I TE ( 1 08, 100)  SYIN, YSC, YRE, HHR, 82 
..RITE(  108,110)  ( PARM  ( I  ) ,  I  « 1, 7  ) 

L9  20  J«l,5 
A  N  *  0  *  0 

A«0 


A -5 


5 


FRR«FR7 

N  »I<  +  1 

FRQUWFRR 
^■TRtiP  1  *FRR 
AXR»*/CZERtj 
i F  ( J-t )  6/7/R 

6  ST6P  6 

7  CALL  I MP ( HD/ S ! GO ) 
fcl«ZZEK/H 

1.J1  U)*Dl 
bGD ( A ) «  S I G 
L6  Tr)  ^ 

8  Li«ani(K> 

SIG«SGu(K) 

9  AXsXZ ( J ) *A*R 
RR»HHR*AKR 

r-i I  sHHl  »AxR 

X2«XX«XX 

bY«Sr I\*AKR 

S  Y  «o  Y ♦HR 

Z  »HR *HR+h I *H I 

P5G«0«C 

L'9  18  I  »l/b 

YY*YZ( I )*YSC#AKR+hR+YRF*AKR  - 

T 1 *YY-SY 

72«Tl*Tl 

RR«T2*x2 

R«5uiRT  (  RR  ) 

T1«YY*SY-HR*2.C 

T2«T 1 *T 1 

RRl2«T?+x2 

*Rl«SuRT(RRl?) 

uZR»T1/RR1 

uZR2«oZR*GZR 

FN«0.5*( 1.0-GZRP) 

CALL  Cvl V2 

EN1»Fnc*vDD1-GZ*VD1 

t\2*FNC*VDD2-uZ*VD2 

R2*Rl  +  (2.0/0.0)  #H*(Gl*Slr.-GZ) 

7EMP»((.;.0/1»0)*(R2”R7) 

TEMP«CEXP(TL.vP) 

7EKP«TF.MP*RZ/R1 
Tt.vlP2«Rl/  H 

T^*TLivlP#(V2”(C»0*  1  •  0  )  *EN?/TE -1P2  > 
7EMP»(D.0/l.O*(  (F\12-RZ2)/(R1+RZ)  ) 
IcitfP-CExPtTEPP) 

TEiiP*TFMP»RZ/Rl 
Tt»P2»ENl/Rl 
TEMP2*TErtP2*(0»0# 1»0> 

T  EMP2* V 1-TEMP2 
TtMP«TtMP»TEMP£ 

T  >%  *  T  «♦(  1 .0/0.0 
7  .vsTa  +  TEMP 
AaC*7wR#TrtR+TkI#TrtI 

18  PS'j«PSG+ABC 

CBP ( K ) »5»0/PSS 
AN»A,\i+1 .0 

FhR»FR7 #10.0** ( AN/ANN) 

IF  (K.LT.301)  GO  T9  5 
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CALL  PLdTIT  (FR3*8eP*301* IPI < J) ) 

CALL  T8v ( J) 

20  CONTINUE 

A ( 1 ) *AkD( 1 ) 
ca  22  I-2*3C 

A-  (  1*1  )  *10  +  1 

22  X ( 1 ) -AkD ( K ) 

X ( 31 ) -AKD ( 301 > 

09  30  u-l*5 

CALL  PlQTIT  (X*Y( 1* J)*31* IP2< J) ) 

3C  CONTINUE 

*RITt  ( 108*31 ) 

31  FORMAT  ( ///16X* 45H9NE-TH1 RD -OCTANE  RAND  EXCESS  ATTENUATION*  DR./ 
lllX*55(lM-)//tlx*4HEREC*18X*?3HRANlGF  DISTANCES  -  FEET./1?X*?HHZ*9X 
2*  3H25C*  7X*  5H  500.  *  5X*  5Hlo00<'  *5X*  5H20C0.  *  5X*  5H4000./1  IX*  4  ( 1H-  )  *  7x*  4 
35 ( 1H« ) / ) 

09  4.1  1-1*31 
C J  32  J»l*6 

32  C.l  (  J )  » 1  0*  0*AL8G1  C  (  Y  (  I  *  J )  ) 

*RITE  (108*34)  X ( I ) * D8 ( 5 ) * DB ( 4 ) * DB ( 3 ) # DB ( 2 ) * DB ( 1 ) 

34  FOkMaT  (bX*FlO»?*5X*5(F8.2*2X)  ) 

4C  CONTINUE 


..kite 

(103*50) 

FORMAT 

( 1  H 1  ) 

u3  TO 

2 

P  CJRMAT 

( 1CX*5HSY 

*£12*5/ 

1 

10x*5hYSC 

*E12.5/ 

A 

1QX*  5HYRE 

*  £  1 2  »  5  /  •' 

a 

10X*5HH 

*F1?»5/  . 

c 

10X*  5hb2 

*E1P.5*3H  *  *E1?.5) 

f  Ufi.M  AT 

( 10X*5HSi< 

*E12.5/ 

1 

1 Oa*  5hY 

*E12.5/ 

2 

1 OA*  5HR 1 

*E12*5/ 

3 

1  OX*  5HRM 

*E12.5/ 

4 

1  Oa*  5hu, 

*E12.5/ 

5 

1  Ox*  5HCK 

*E12.5/ 

6 

1  OX*  5nRZ 

*  E 1 2  •  5  ) 

999  STOP  999 
END 


subroutine  tovu) 

COMMON  SPARE (1204) 

COMMON  XI (301 )*Y1(301 ) * X? ( 31 )*Y2(31*5) 

SUM-l.o/YK 1 )+l.C/Yl(2)+1.0/Yl ( 3 ) +1 . 0/ Yl £  4 ) +1 . 0/Yl (5) +1.0/Y1 ( 6 > 
Y2(1*J)«6. O/SUM 

SUM- 1 • C/Yl ( 300 ) +1 . O/Yl ( 299 )+l » /Y1 <298>+l*/Yl ( 297 ) +1 • /V 1 ( 296 ) +1 • 0/Y 
11(301) 

Y2(31*J)-6. O/SUM 
CO  10  1-2*30 
SUM-0.0 
KE-U10-5 
KS«KE"8 
09  6  <«KS*  KE 
6  SUM-SUh+l.O/YKK) 

SUM-SUM  +  2.Q/(Y1 (KS-1 )+Yl (KE  +  1  )  i 
1C  Y2 ( I  *  J ) « 10.0/SUM 
RETURN 
END 
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subroutine  cviva 
COMPLEX  TEBPA*  TEMPB/ H 

CBl'PcEx  u2/Gl>bl/B2/SIG/SI2/Vt/VCl/VDDl/V2/VD2/VDD2/G72/Gl?/TEKP 
CaMM9N/$vv/G2»8t/B2/Vl/ VC/l/VCDl/V2/VD2/VDD2/Gl/HjSIG/TW6PI 
C I y£NS I 6\  SO ( 2 ) 

tjUi  VALENCE  (SC/SIG/Sift)/  ( SO ( 2 ) / S 1 1  )/  (G2R/G7) 

c  *** 

aS* AbS ( S I  I ) 

IF  (XS-1.0E-9)  ICO/ 100' lu8 
ICO  x5«ArS ( 3  IP i 

AS«Ar>S  (  XS»  1  •  0  ) 

IF  (xS-l.OE-9)  102/ 102/ 108 
1C2  vL)l»(0. 0/0.0) 

v  UD 1 «  v  C 1 

V 1 «  C  ( 1.0/0»0)-Bl )/( (1*0/ 0*0) *81 ) 

V2« (4.0/C*0)*Bl*(G2-e2)/<G2+b2)/( (1.0/0*C)+Rl )/< ( l*0/0*0)+Bl ) 
xS*AdS <  GZR ) 

I F ( X3-2  • 0E* 1 )  10A/10A/106 
1C*  Vu2»VG1 
vOD2» VC2 
f-ET  URN 

1C6  VU2»(8*0/0*G )*el*B2/( ( 1 *0/ 0*C ) +B1 ) / ( ( 1 » 0/ 0*0 > *Bl ) / ( GZ+B2 ) / ( GZ+B? ) 

V002« VC2# ( «2*  0/ 0*  0 ) / (GZ+B2 ) 
ulso2 
oO  T9  ?00 
1L8  CSNTINuE 

SI2*3IG*SIG 

o22«G2*GZ 

oUSI2*G;2-(1.0/C.0) 
ulsCGQkT (SI) 
ul«Gl/SlG 
C1?«G1*G1 

vlx (o2-Bl#Gl ) / ( G2+9 1 *G1 ) 

\,01*(2.0/0.0)*bl*(G12*SI2-G22)/SI2/Gl/{u2+Bl*Gl  )/(GZ+Bl*G1  ) 
vU?  =  -vrQl*(o22+(?.0/0»0)*Gl2*S12+<3»n/0«0)#G7*Gl»Bl  ) / ( S I2*G12*  (  G7+B 
11*31 )  ) 
v0Dl«V02 
TEMP* ( uZ+b 1* G1 ) 

TEMP«TEMP#TEf^P 

v  2*  (  31  -32  )  »  ( * .  0/  C*  0  )  *81  *GZ*G1  /  (  TE*lP*  ( G1+S2  )  ) 

X 3  =  A  iS  (  GZ R  ) 

IF  US-1.0E-9)  lfcO/ 160/ 1 7C 
ItO  vu?*(0* 0/0*0) 
v  Jl)2*Vu2 

^cTjRN 

17C  V J2* ( 8 • 0/ 0*0 ) *B1 *B2*G22/ ( SI 2*TEMP* ( G1+B2 ) # ( Gl+62 ) ) 
v  J2*  ,/D2  +  V2*Vl*( ( 1*0/ 0*0 )-SI2)/Sl2/Gl 2/62 

v JD2»( ( 1*0/0.0)-SI2)/SI2/G12/GZ*( ( VI *vC2+ VD 1  * V2 ) - V1*V2* { S 1 2*Gl2+3Z 
12* (2*0/ 0*0) )/SI2/GZ/Gl2) 

TEmP»GZ/( ( o2*Bl*Gl ) *S I G* ( G1 +82  )  ) 

7tMP»TEMP*TEMP 

v J02*VuD2+( 16*0/0*0) *B1 *B2*TEMP* ( ( 1 .0/ 0* 0 j /GZ- ( SI 2*G1+B1 *GZ ) /Si 2/G 
ll/(G2+bl*Gl J-G2/5I2/G1/ (G1+B2)  ) 

20C  TEriPA»G2*(0*0/2*0)/SIG/Gl 

TdpPa»TEPPA*TEKPA*h+(O.G/  2* 0 ) * ( S 1 2- ( 1 . 0/  C. C ) ) /G12/G1 /S 1 2/SI G 

VDD2»VCD2/H 

V02» VD2/H 

v0D2»VD02+(2*0/0.0)*TEMPA#vD2+TEMPB#V? 

V02« VD2+V2*T£MPA 

RETURN 

END 
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SUBROUTINE  PL0TIT  (X#Y#N#IP) 

DIMENSION  X(1)#Y(1) 

C0MM8N/PLUT /XST #  XAX IS# XDECADE#  YAXlS# YDECADE# I Q#  YMIN 
C0MM0N/AX I S/TT ( 8 ) 

DATA  TT/'u*  0458#  0*  0369#  0*1 549# 0*2218#  0*301  #0*3979#  0*5229#  0.699/ 

DATA  XST#  XDECADE#  YDECADE# I Q/3.0# 2*0# 1 .5# 5/ 

DATA  YAXlS#YMIN/7.5#-10*0/ 

YSftYDECADE/10.0 
IF  (IP-1>  230# 10# 17Q 

10  IF  (ABS(X(l)*l*0)-0. 00001)  20# 20# 30 
20  XSCALE-I.O 
60  T0  60 

30  XM«AL0G1O(X( 1 )  ) 

IF  (X(l>*1.0>  40# 20# 50 
40  I»IFIX(XM)-1 

XSCALE-1. 0/10. 0**1 
G0  T0  60 

50  I « IF  IX ( XM+0.5 ) 

XSCALE«1. 0/10. 0**1 
6L  CALL  P10T  ( XST #  0*0#  25 ) 

I«1FIX(AL0G1O(X(N)*XSCALE)+1.O) 

XAXlS*FL8AT< I )*XDECADE 

XST  «XAX IS  +  3‘0 

CALL  PLOT  (XAXlS#0*0#2) 

CALL  PL0T  ( XAXl S#  YAX I S#  2 ) 

CALL  PL0T  (0.0# YAXIS#2) 

CALL  PL0T  ( 0. 0# 0*0# 2 ) 

XX«X0ECADE 

L0  80  K  »  1  #  I 

L8  7o  J«l#8 

X0*XX-TT ( 9" J ) *XDECADE 

CALL  PL0T  ( XD# 0. 0625# 1 ) 

70  CALL  PLOT  (XD#0.G#2) 

CALL  PL0T  (XX#0. 125# 1 ) 

CALL  PLOT  ( XX#  C»  C#  2 ) 

8C  XX»XX+XDECADE 
XX»X AX  I S 
TYsC.O 

CO  110  K«l#  IQ 
YD«YY 

LO  90  J»l#4 
YD« YD+YS 

CALL  PL0T  (XX# YD#  1  ) 

9C  CALL  PLOT  ( XX-O. 050# YD# 2 ) 

YD«YO+YS 

CALL  PLOT  ( XX#  YD# 1 ) 

CALL  PLOT  (XX-O. 1# YD#2) 

L 0  100  J* 1# 4 
YD«YD+YS 

CALL  PLOT  (XX#YD#i) 

ICO  CALL  PLOT  ( XX-O. 050# YD# 2 ) 

YY.YY+YDECADE 
CALL  PL0T  (XX#YY#1) 

CALL  PLOT  ( XX-O. 150#  YY#  2 ) 

110  C0NTINUE 
YY«YAXIS 
08  130  K» 1#  I 
DO  120  J* 1#  8 
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XCaXX-TT<J)«XDECADE 
CALL  Pl8T (AD4YY4 1 ) 

CALL  PL8T(XD,YY-C«0625*2> 

XXsXa-xDECADE 

CALL  Pl6T(xX,YY,1) 

CALL  PL8T(XX>YY-C«125#2) 

D8  160  K«1/IC 
YD«YY 

L8  140  Ja 1 4 4 
YOaYD-YS 

CALL  PL8T  (O.O/YD/l) 

CALL  PL8T  (0*0504 YD/ 2 ) 

YJaYO-YS 

CALL  PL8T  (O.O/YD/l) 

CALL  PL8T  (0*1/  YD#  2 ) 

C0  150  J» 1 4  4 
YOaYD-YS 

CALL  Pl6T  (O.O4YD4I) 

CALL  PL0T  ( 0*050/ YD4 2 ) 

YYaYY-YDECACE 

CALL  PL8T  (O.O4YY4I} 

CALL  PL0T  ( 0* 1504  YY4  2 ) 

CONTINUE 

XX«AL0Ol 0 ( X ( 1 )*XSCALE)*XCECADE 

YYalj«o*AL9GlO(Y( 1 ) l-YMIN 

IF  ( Y  Y )  180/1904190 

YY.G.O 

YYaYY*YS 

IF  <YY.GT*9*0)  YYa9.0 
CALL  PL0T  <XX4 YY4 1 ) 

Ld  220  I«24N 

XXa AlBoIO ( X ( I )*XSCALE)*XDECADF 
YY«iy.C*AL8G10(  YU  )  >-YMl\ 

IF  </Y)  20042104210 
YYaG.O 
Y  Y*YY*YS 

IF  (YY.GT«5*G)  YY«9*0 
CALL  PL6T  (XX4YY42) 

C0NTIi\uE 

PETUSN 

END 


A-10 


I 


I 

I 


J 

\ 


SUBROUTINE  IMP(ACDaSIGD) 

C 

CQMPuEa  BC#  XC#  ZC 
C 

COMMON/LIMI TS/CZER9 
r 

common  spareii204) 

COMMON  FREQ (301 )* DUMMY <4S2)jSK# Y,R1,RMaQ/CK* R?*Ut DUO 
COMMON  ZCD(2)*BCD(2)/ XCD ( 2 ) 

C  *# 

DIMENSION  ACD(2)iSIGD(2) 

C 

EQUIVALENCE  (BCDxBCiBR)#  <BCD(?)iBI ), (XCD/XCaXR)> (XCD(2)*XI ) 
EQUIVALENCE  (ZCD<7C>ZR)j (ZCD(?)#Zl  ) 

c  ** 

IF  ( Ck-G#2U« 0  )  100>100>1?G 
ICO  aK»1«C 

a  I *-Rl  /RZ/SK/l* 
xC*CSGRT(XC) 
cRsO • 0 

bI*RZ*SK*Y/C</1.4 
cl «w*SQRT (BI  ) 

LC*BC*XC 
110  zRsO.O 

Zl=*CK/w/Y*l#4 

ZC*ZC*BC 

CMeW/BI 

S I GD ( 1 )«CZER8/CM 
SIG0(2)«BR«CZER8/W 
CO  TO  130 

120  TERM1«1.Q+RZ*(Sk-1«0)/RM 

TEkM12»TERM1 *TERM1 
T£RM2«RM*W 
TERM22«TFRM2#TERM2 
A 12*Rl*Rl 

7tRM3»TERM22*TERMl2 
1 3T4*1  • 0+R 12/TERM3 
kZk«RZ#Sx 

TOP* 1 »0+Rl2* ( Y+RM/RZK  J /TERM3 
f-0T2*l.C+RZ*(SK-l»C)  /RM 
R  31 *RZ«*T0P/B6T4 
KRB1*1.0+RZ*(SK-1.0)/RM 
R  Ri-  lsRRbl*RRBl 
c0T4*6OT  4#RR61 

A  Rt 1 *Rl * ( 1  • 0»RZ* ( 1.0-Y I/RMJ/R0T4 
A  K  *  R  b  1 
X  I *-RRbl/w 
xC*CSGRT ( XC ) 
r;I*Y/CK/l  .4 
b I *30RT ( B I  ) 
c I *b I  * a 
E.  K  ■  0  •  0 
t,C«BC*xC 
CO  TO  110 
130  Continue 
AGO ( 1 ) «ZR 
ACD ( 2 ) *Z I 
RETURN 
END 
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